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1.  IKTRODUCTION 


The  research  summarized  in  this  report  is  a  continuation  of  that 
reported  earlier  {Ref.  1-11).  The  primary  objective  is  to  understand 
and  learn  how  to  control  combustion  instability  in  solid  rocket  motors, 
especially  in  the  low  frequency  range  (up  to  about  500  cps).  This  ob¬ 
jective  was  stimulated  by  the  Increasing  interest  in  very  large  booster 
motors  and  low  L*  space  motors.  Laboratory  tests  and  theoretical  con¬ 
siderations  indicated  a  significant  risk  of  low  frequency  instability, 
and  the  impact  of  such  problems  on  development  programs  promised  to  be 
very  costly  unless  better  laboratory  tests  and  understanding  of  the 
phenomenon  were  available. 

Preliminary  studies  were  concerned  with  development  of  several 
laboratory  techniques  for  the  study  of  low  frequency  instability  (LFI), 
and  use  of  these  methods  to  screen  the  effect  of  propellant  variables 
on  LFI  (Ref.  1-4).  This  work  provided  a  preliminary  basis  for  assess¬ 
ment  of  the  LFI  problem  in  rocket  motors,  and  provided  experimental 
methods  that  could  be  applied  when  problems  arose  in  development  pro¬ 
grams.  However,  it  became  evident  that  an  adequate  approach  to  the  LFI 
problem  was  dependent  on  a  better  understanding  of  the  detailed  aspects 
of  propellant  combustion,  and  several  studies  (Ref.  4-7)  were  initiated 
which  were  of  a  more  fundamental  nature,  related  to  ingredient  decom¬ 
position,  combustion  zone  structure,  and  acoustic  attenuation  processes 
(Ref  8-10).  These  studies  were  continued  along  with  the  LFI  studies. 

The  fundamental  studies  inc  uded  decomposition  of  ammonium  per¬ 
chlorate  (AP)  by  differential  thermal  analysis  (DTA)  and  by  observation 
of  nucleatlon  and  spread  of  decomposition  sites  on  large,  single,  high 
purity  crystals;  observation  of  AP  deflagration  by  high  speed  photogra¬ 
phy  and  study  of  quenched  samples  of  large  single  crystals;  study  of 
combustion  zone  structure  by  photography  and  quenching  of  "sandwiches” 
of  AP  and  binder;  photographic  observation  of  accumulation,  agglomer¬ 
ation  and  combustion  of  powdered  metal  fuels  in  propellant  cumbustion; 
hot  plate  heating  and  agglomeration  mechanism  of  metals,  snd  cold  flow 
studies  of  acoustic  attenuation  ip.  model  combustors.  These  areas  of 
research  have  been  presented  extensively  in  previous  reports  and  publi¬ 
cations.  During  the  past  year  it  has  been  necesssry  to  drop  some  of 
these  areas  of  investigation  due  to  reductions  in  funding.  Specifi¬ 
cally,  work  in  the  following  areas  was  dropped  or  transferred  to  other 
programs:  AP  decomposition,  "sandwich  burning",  and  experimental  work 
on  acoustic  attenuation.  The  areas  of  continuing  effort  are  indicated 
by  the  section  headings  in  the  report. 
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NOMENCLATURE 
E  T 

A  Nondlmensional  parameter,  (l  -  — ) 

s  s 

2/3 

8  Parameter  in  the  granular  diffusion  flame  correlation  p/r  “  a+bp 

E  T 

B  Nondlmensional  parameter  - 1/2)  (l  -  ) 

s  s 

2/3 

b  Parameter  in  the  granular  diffusion  flame  correlation  p/r  -  a+bp 
c*  Effective  discharge  velocity  of  nozzle 
Eg  Activation  energy  of  solid 
L*  Volume  of  chamber/area  of  nozzle  throat 
m  Mass  burning  rate 

n  Pressure  exponent  from  r  »  Cpn 

p  Pressure 

R  Magnitude  of  combustion  response  to  pressure  disturbance, 
ft.  -  (m/m) /(p/p) 

r  Burning  rate 

o  Growth  rate  of  pressure  oscillstior.s 

P  "Constant”  value  between  1  and  (r+l)/2r,  depending  on  extent  of 
temperature  oscillation  in  combustor 

k  Thermal  diffusivity 

X  Characteristic  root  of  solid  phase  energy 

p/e  Ratic  of  mass  perturbation  to  pressure  perturbation.  Also  desig¬ 
nated  the  propellant  response  function 

T  Time  lead  of  combustion  oscillation  relative  to  perturbing  pres¬ 
sure  oscillation 
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Much  of  the  work  on  the  present  program  has  seen  publication  In  two 
or  more  forms,  e.g.,  as  progress  reports,  sunmary  reports  of  particular 
areas  of  work,  and  as  open  literature  publications.  In  the  Interests  of 
economy,  it  was  decided  that  the  present  report  would  present  only  brief 
resumes  of  those  areas  In  which  complete  publication  In  some  other  form 
was  available  on  a  comparable  time  scale. 
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2.  AMMONIUM  PERCHLORATE  DEFLAGRATION 


Wien  attention  was  directed  to  the  structure  of  the  combustion  zone, 
it  seemed  appropriate  to  concentrate  concurrently  on  the  deflagration 
behavior  of  AP,  since  it  is  usually  t!.“  dominant  propellant  ingredient. 

On  burning  and  quenching  single  crystals  of  AP,  it  was  quickly  found 
that  prevailing  views  regarding  the  deflagration  were  extremely  naive. 

The  burning  surface  exhibited  certain  reproducible  irregularities,  that 
were  dependent  on  pressure  (Ref.  4,  7,  12,  13  and  14).  A  previously  un¬ 
reported  liquid  layer  was  observed  (Ref.  4,  5,  7,  12,  13  and  14),  which 
appeared  to  encompass  a  significant  portion  of  the  rate-controlling 
reactions.  Although  such  a  layer  is  not  embodied  in  any  analytical 
models,  its  role  seems  to  be  crucial,  and  its  presence  resolves  aome 
earlier  controversies  regarding  the  rates  and  relative  importance  of 
solid,  surface  and  gas  phase  reactions.  During  the  present  reporting 
period,  these  studies  have  been  extended  to  higher  pressures,  and  have 
revealed  even  more  complicated  behavior  (Ref.  14).  The  investigations 
of  the  deflagration  of  single  crystals  were  carried  out  to  ascertain; 

(1)  the  variation  of  the  deflagration  rate  as  a  function  of  pressure, 

(2)  the  surface  structure  of  the  deflagrating  samples,  and  (3)  the  sub¬ 
surface  profile  associated  with  the  deflagration.  From  theae  observa¬ 
tions  energy  transfer  modes  were  postulated. 

Data  concerning  the  variation  of  deflagration  rate  versus  pressure 
from  our  investigations  and  those  of  other  laboratories  are  compared  in 
Fig.  2.1.  Our  results,  obtained  using  cinemicrophotography ,  confirm  the 
trend  described  by  Bobolev,  et  al,  and  Glaskova.  Come  differences  be¬ 
tween  their  data  and  ours  might  be  expected  aa  they  used  pressed  pellets 
and  we  used  single  crystals  of  AP. 

We  observed  that  the  single  crystal  data  define  a  curve  which  nat¬ 
urally  appears  to  be  divided  into  four  regimes  (Fig.  2.2).  Not  only  do 
these  divisions  apply  to  the  deflagration  rate  curve  but  also  to  the 
results  of  the  surface  and  aubaurface  structures  of  the  deflagrating 
samples.  The  behavior  la  summarized  in  Table  2.1.  The  four  regimes 
characterize  the  deflagration  of  AP  crystals  at  pressures  from  300  to 
10,000  psla,  thua  permitting  predictions  of  AP  behavior  during  defla¬ 
gration  from  a  knowledge  of  the  pressure. 

Scanning  electron  microscope  (SEM)  micrographs  taken  of  samples 
from  Regime  1  Indicate  thst  the  surface  was  covered  with  a  froth:  gas 
evolved  through  a  liquid  layer  (Fig.  2.3).  This  obaervation  permitted 
a  choice  to  be  made  between  two  competing  hypotheses.  The  widely 
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FIG,  2.1.  Deflagration  Rate  of  Amronium 
Perchlorate  in  Nitrogen  Atmosphere. 
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FIG.  2.2.  Deflagration  Rate  of  Ammonium  Per¬ 
chlorate  Single  Crystals  in  Nitrogen  Atmosphere 


FIG.  2.3.  Scanning  Electron  Micrographs  of  Ammonium 
Perchlorate  Single  Crystals  Which  Were  Quenched  While 
Undergoing  Deflagration  at  the  Pressures  of  Regime  I. 


accepted  assumption  that  the  mechanism  of  AP  deflagration  was  the  same 
as  decomposition  {i.e.,  sublimation)  waa  not  seriously  questioned  until 
Hightower  (Ref.  4,  5,  and  12)  observed  a  cooled  froth  on  the  surface  of 
quenched  samples  which  he  examined  using  an  optical  microscope.  Scan¬ 
ning  electron  microscope  analysis  revealed  the  details  of  thia  froth 
and  it  was  apparent  that  a  melt  had  been  present  on  the  surface  of  the 
deflagrating  sample.  Such  observations  proved  that  the  sublimation 
hypothesis  was  inaccurate  to  describe  the  deflagration.  The  froth  ou 
the  AP  crystal  surface  was  shown  to  be  roughly  1  to  5  microns  thick  on 
the  300  and  600  psia  samples  and  present  only  in  patches  cn  samples 
burned  at  800-900  psia.  Examination  of  the  profile  showed  that  this 
froth,  the  thickness  of  which  decreased  with  increased  pressure, 
covered  a  layer  of  material  which  had  been  cubic  phase  AP,  which  In 
turn  covered  tho  original  orthorhombic  phase.  The  well  defined  inter¬ 
face  between  the  cubic  and  orthorhombic  phases  indicates  the  location 
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of  the  242°C  isotherm  within  the  crystal.  Knowledge  of  the  position  of 
this  Isotherm  and  the  regression  rate  as  a  function  of  pressure  allows 
a  calculation  of  the  surface  temperature  of  deflagrating  crystal 
(ief .  15) .  Results  obtained  by  using  the  SEM  data  were  more  precise 
than  those  Beckstead  and  Hightower  obtained  using  the  optical  micro¬ 
scope,  and  were  in  agreement  with  the  "best"  line  drawn  through  their 
data. 

Samples  from  Regime  II  differed  from  those  of  Regime  I  in  that  no 
longer  was  there  a  liquid  on  the  surface,  rather  the  surface  was  charac¬ 
terized  by  a  pa 'tern  of  ridf.es  and  valleys  (Fig.  2.4a)  with  sites  of 
increased  activity  at  the  bottom  of  the  valleys  (Fig.  2.4b).  Ridge 
length  was  found  to  vary  inversely  with  the  pressure  (over  500  microns 
long  at  1000  psia  and  less  than  100  microns  at  1800  psia) .  Jest  as  the 
ridge  structure  was  pressure  depandent  so  were  the  sites  of  increased 
activity  (reaction)  at  the  bottom  of  the  valleys.  The  samples  quenched 
at  1000  psia  show  frothy  residue  in  the  valleys.  The  samples  quenched 
at  1200  psia  show  frothy  material,  but  also  show  holes  of  complex  geome¬ 
try  indicative  of  localized  activity.  At  1500  psia  an  array  of  needle- 
like  structures  appeared  In  some  of  the  holes.  This  feature  becomes 
Increasingly  evident  with  increasing  pressure. 

The  profiles  of  F.egime  II  samples  present  evidence  of  the  complex 
nature  of  the  combustion.  The  cubic  phase  thickness  was  greater  under 
the  ridges  than  under  the  valleys  (Fig.  2.4c).  This  difference  in 
thickness  indicates  that  the  regression  rate  of  the  ridges  was  less  than 
that  of  the  valleys,  that  the  temperature  gradient  through  the  ridges 
was  not  as  steep  as  under  the  valleys,  and  that  heat  was  transferred  in 
three  diversions  (including  conduction  of  heat  in  direction  perpendicu¬ 
lar  to  the  direction  of  the  average  surface  regression) .  This  inter¬ 
pretation  would  explain  why  the  ridges  frequently  peeled  off  from  the 
surface  without  burning,  as  seen  in  the  high  speed  motion  pictures. 

Figure  2.4d  indicates  that  the  activity  sites  of  the  valley  extended 
below  the  surface  for  a  depth  almost  equivalent  to  the  cubic  phase  . 
thickness. 

The  surfaces  of  the  Regime  III  samples  show  pronounced  depressions 
(Fig.  2.5)  which  have  extensive  areas  of  needles.  The  micrographs  of 
the  profiles  indicate  that  the  needles  extend  200  to  300  microns  from 
the  surface  of  the  pocket.  Studies  using  cinemicrophotography  identi¬ 
fied  the  origin  of  these  depressions  and  helped  explain  why  the  needles 
exist  only  in  the  depressions.  These  high  apeed  motion  pictures  show 
localized  orange  flamelets,  which  change  location  as  the  burn  progresses, 
over  the  surface.  Parts  of  the  surface  were  seen  to  regress  locally 
(cat.1  ’  the  depressions  seen  in  the  micrographs)  while  adjoining  areas 
were  temporarily  stationary,  with  the  flame  standing  over  the  recessed 
areas.  Because  the  flame  was  observed  to  be  in  the  depression,  the  lip 
seen  in  Fig.  2.5  provides  an  opportunity  to  view  a  portion  of  the  sample 
which  was  partially  exposed  to  the  flame  and  at  the  same  time  view  a 
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FIG.  2.4.  Scanning  Electron  Micrographs  of  Ammonium 
Ferchlorate  Single  Crystals  Which  Were  Quenched  While 
Undergoing  Deflagration  st  the  Pressures  of  Regime  IX. 


portion  which  wss  shielded  from  the  flame.  Xhst  portion  which  tha  flame 
contacted  exhibits  the  needle  structure,  while  the  shielded  portion  ex¬ 
hibits  s  different  surface  structure.  Thus  by  combining  the  SEM  re¬ 
sults  with  those  obtained  using  cinemlcrophotography  it  is  possible  to 
show  that  the  needle  structures  sre  indicative  of  proximity  to  the  flame. 
Scanning  electron  microscope  micrographs  taken  of  Regime  IV  samples  show 
that  needles  covered  the  entire  surface  (Fig.  2.6)  and  the  high  speed 
motion  pictures  show  a  pronounced  orsnge  flame  adjacent  to  the  burning 
surface,  thus  confirming  the  above  hypothesis. 
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From  these  observations  we  hypothesized  ss  to  how  energy  wss  trans¬ 
ferred  to  the  unreseted  AP.  At  the  low  pressures  of  Regime  I  we  have 
seen  that  the  surface  is  covered  with  s  froth.  It  Is  our  hypothesis 
that  decomposition  occurs  either  within  or  beneath  the  melt  and  that 
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FIG.  2.5.  Scanning  Electron  Micrographs  of  Ammonium 
Perchlorate  Single  Crystals  Which  Were  Quenched  While 
Undergoing  Deflagration  at  the  Pressures  of  Regime  III. 


exothermically  resctlng  gss  products  are  entrapped  within  the  melt. 

The  heat  produced  by  the  exothermic  reaction  is  transferred  by  conduc¬ 
tion  through  the  froth  to  the  unreacted  solid,  thereby  providing  the 
energy  required  for  further  decomposition.  Further  exothermic  reaction 
prohsbly  occurs  in  the  gss  above  the  surface,  althcugh  the  heat  trsns- 
fev  from  those  resctlons  may  be  secondary  in  importance.  As  the  pres¬ 
sure  is  increased  the  liquid  becomes  so  thin  that  it  ruptures  into  drop¬ 
lets  and  the  high  reaction  rate  in  these  droplets  is  responsible  for  pro¬ 
ducing  the  valleys .  The  reactions  become  concentrated  within  the  valleys 
and  three-dimensional  heating  occurs  within  the  solid. 

While  exoth:rmic  reactions  within  the  froth  and  sctivity  sites  within 
the  valleys  are  responsible  for  most  of  the  energy  trsnsfer  occurring 
in  Regimes  I  and  II,  energy  is  primarily  transferred  to  Regime  III  and 
IV  samples  by  conduction  from  the  gas  flsme.  In  Regime  III  the  unsteady 
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FIG.  2.6.  Scanning  Electron  Micrographs  of  Ammonium 
Perchlorate  Single  Crystals  Which  Were  Quenched  While 
Undergoing  Deflagration  at  the  Pressures  of  Regime  IV. 


flame  cauaed  an  areawise  "atop  and  go"  regression  of  the  surface.  In 
Regime  IV  the  flame  la  steady  and  the  deflagration  rate  once  again  in- 
creaaea. 

A  more  complete  deacription  of  the  recent  work  under  this  program 
is  contained  in  Ref.  6,  14,  16  and  17.  From  the  results  it  is  clear 
that  AP  deflagration  is  a  very  complex  process,  quite  different  from 
the  assumptions  of  analytical  models,  and  possibly  amenable  to  modifi¬ 
cations  not  previously  suspected.  It  is  interesting  to  note  that  most 
analytical  modela  of  composite  solid  propellant  combustion  rest  on  in¬ 
correct  assumptions  regarding  AP  deflagration,  including  AP  deflagra¬ 
tion  rate  curves  that  completely  Ignore  the  drop  in  burning  rate  above 
2000  psi  Thus  the  propellant  burning  theories  are  gravely  in  error, 
even  if  they  do  (with  adjustment  of  undetermined  parameters)  successfully 
correlate  experimental  results. 
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3.  PROPELLANT  COMBUSTION  MECHANISM 

The  most  conspicuous  problem  in  combustion  instability  research 
has  been  the  continuing  ignorance  of  the  nature  of  the  combustion  zone 
whose  perturbationa  cause  instability.  During  the  present  reporting 
period,  several  investigations  were  directed  towards  better  understand¬ 
ing  of  cteady-atate  combustion  and  combustion  zone  structure. 

3.1.  ONE-DIMENSIONAL  ANALYSIS  OF  COMBUSTION  ZONE  STRUCTURE 

Most  anslyticsl  models  of  steady-state  burning  of  solid  propellants 
can  be  referred  to  as  "complete"  in  a  mathematical  sense,  in  that  solu¬ 
tions  are  entirely  a  consequence  of  the  set  of  equations  in  the  model 
(with  a  few  undetermined  constants).  Such  modela  are  both  complex  and 
naive.  They  are  complex  in  that  even  simplified  representations  of  the 
combuation  zone  behavior  are  analytically  complicated.  They  are  naive 
in  thst  the  simplified  representations  ignore  important  aspects  of  the 
real  process  (such  as  three-dimensionality,  Innomogeneity  of  the  pro¬ 
pellant)  . 

In  attacking  the  problem  of  combustion  zone  structure,  it  appeared 
that  one-dimensional  models  were  really  only  computational  schemes  to 
aid  physical  insight,  and  as  such  should  be  formulated  with  a  maximum 
flexibility  to  accommodate  experimentally  observable  information  such 
aa  combustion  zone  thickness  and  burning  rate-pressure  data.  With  this 
philosophy  in  mind,  a  one-dimensional  model  was  constructed  which  was 
not  "complete",  and  could  accommodate  input  data  on  combustion  zone 
locstion  and  thickness,  burning  rate,  etc. 

The  model  Is  an  elementary  one  in  which  an  assumed  surface  pyroly- 
ais  is  followed  by  a  gas  phase  reaction.  The  gas  phase  reaction  rate 
la  assumed  to  be  constant  in  the  reaction  zone,  which  starts  at  some 
parametrically  specified  distance  from  the  surface.  The  solution  used 
here  is  aort  of  an  Inverse  one  in  which  all  required  thermodynamic, 
kinetic  and  transport  properties  are  specified  as  well  as  the  burning 
rate-presaure  curve. 

Calculations  based  on  observed  burning  rate-pressure  functions 
yiald  values  for  the  heat  release,  flame  thickness  and  stand-off  dis¬ 
tance  required  to  give  the  specified  burning  rate.  Using  plausible 
values  for  thermodynamic,  kinetic  and  transport  properties,  the  results 
for  flame  thickness,  etc.,  are  well  within  the  ranges  observed.  It 
appears  then  that  this  approach  may  be  useful  for  examining  the  gross 
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FIG.  3.1.  Burning  Rate  Curves  for  Polyurethane' 
Ammonium  Perchlorate  Propellants . 
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effects  of  some  of  the  parameters,  auch  aa  surface  activation  energy, 
and  thereby  helping  relate  burning  rates  to  cnemical  composition.  The 
simplicity  of  the  analysis  enables  one  to  perceive  easily  the  physical 
origin  of  the  results,  a  feature  which  is  often  not  true  In  more  sophis¬ 
ticated  models.  The  results  will  be  useful  in  atudlea  of  more  complex 
models,  be  of  value  for  examining  limiting  conditions  for  time  dependent 
models,  and  be  helpful  in  design  of  experiments  to  measure  combustion 
zone  structure.  Further,  the  calculations  will  help  to  recognize  the 
conditions  under  which  one-dimenslonal  models  are  useful,  and  conditions 
under  which  tbree-dimenslonal  modeling  is  important. 


Results  to  date  are  summarized  in  detail  in  Ref.  22  and  23. 
3.2.  BURNING  RATE  OF  "B1M0DAL"  PROPELLANTS 


It  was  observed  in  bulk  mode  instability  tests  that  the  trend  of 
Instability  behavior  was  different  for  propellants  made  with  bimodal 
oxidizer  particle  size  distribution  than  for  propellants  with  unimodal 
particle  sizes.  In  correlating  the  test  results,  it  Is  customary  to 
use  the  steady-state  burning  rate  as  a  normalizing  variable  In  the  fre¬ 
quency  scale,  and  burning  rate  determinations  were  made  accordingly. 

It  was  observed  that  the  dependence  of  burning  rate  on  pressure  of  uni¬ 
modal  propellants  was  quite  different  from  the  r  (p)  function  for  bi¬ 
modal  propellants.  In  hopes  that  understanding  of  the  contrasting 
steady-state  behavior  might  provide  some  clue  to  the  contrasting  per¬ 
turbation  behavior,  a  careful  study  of  the  burning  rate  was  made  of 
certain  unimodal  and  bimodal  formulations  used  In  L*  burner  testa. 

The  strand  burning  rates  of  two  classes  of  propellanta  were  de¬ 
termined  as  a  function  of  pressure.  One  series  of  propellanta  consisted 
of  polyurethane  and  bimodal  A1  (PU-AP)  and  the  other  consisted  of 
carboxyl-terminated  polybutadiene  and  bimodal  AP  (CTPB-AP).  The  formu¬ 
lations  for  the  propellants  are  given  In  Table  3.1. 

The  data  for  the  PU-AP  propellants  are  presented  In  Fig.  3.1.  The 
trend  of  higher  burning  rates  for  propellants  having  smaller  oxidizer 
particles  Is  spparent  (compare  the  rates  for  A-146  versus  thoae  for 
A-149).  Also  it  la  seen  that  many  of  these  propellants  extinguish  at 
high  pressure:  the  approximate  upper  deflagration  limit  for  1/4-  by 
1/4-inch  atrands  denoted  by  an  x  in  the  figure.  The  mechanism  responsi¬ 
ble  for  this  extinguishment  is  discussed  in  Section  3.3.  The  meaa- 
burning  rate  and  extinguishment  characteristics  of  the  PU-AP  can  be 
altered  by  Inclusion  of  as  little  ar  IX  of  a  catalyst  such  as  n-butyl 
ferrocene  (Fig.  3.2). 

The  burning  rate  curvea  for  the  CTPB-AP  propellants  are  presented 
In  Fig.  3.3.  These  propellants  do  not  exhibit  the  mesa  characteristic 
nor  the  extinguishment  seen  for  the  PU-AP  propellants.  The  SEM  micro¬ 
graphs  presented  in  Section  3.3  have  shown  that  while  the  polyurethane 
binder  melts,  the  CTPB  binder  gasifies  without  melting.  It  aeems  likely 
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TABLE  3.1.  Composition  of  Research  Propellants 


Ingredients  and  weight  percent 


Propellant 

designation 

Ammonium  Perchlorate 

Binder 

Other 

A-146 

37.5? 

15u 

37.5?  30u 

25?  polyurethane 

.... 

A-148 

37.0? 

15u 

37.0S  200u 

25?  polyurethane 

1?  carbon  black. 

A»lb9 

37.0? 

90u 

37.0?  600p 

25?  polyurethane 

IS  carbon  black 

A-151 

37.0S 

45u 

37.0?  200p 

25?  polyurethane 

i  IS  carbon  black 

A-155 

37.0S 

l*5u 

37.0?  400jj 

25?  polyurethane 

1?  carbon  black 

JU156 

** 

CO 

H 

U"\ 

15v 

22.2?  200p 

25?  polyurethane 

1?  carbon  black 

A-15T 

22.2S 

15u 

51.6?  20 0u 

25?  polyurethane 

1?  carbon  black 

A-158 

37.2S 

90u 

37.5?  6o0p 

25?  polyurethane 

(•«  t 

A-159 

36.0S 

15u 

36.0?  200p 

25?  polyurethane 

1?  carbon  black 
1?  n-butyl 
ferrocene 

A- 160 

36.0? 

l*5u 

36.0?  bOOu 

25?  polyurethane 

IS  carbon  black 
IS  n-butyl 
ferrocene 

A-I67 

37.0? 

15u 

37.0?  80U 

25?  CTPB 

1?  carbon  black 

A-168 

37.0? 

l5u 

37.0S  200 p 

25?  CTPB 

carbon  black 

A- 169 

37.0? 

90u 

37.0?  60  Op 

25?  CTPB 

1?  carbon  black 

A- 170 

37.0S 

i*5u 

37.0?  200p 

25?  CTPB 

1?  carbon  black 

A- 171 

37.0? 

l*5u 

37.0?  kOOp 

25?  CTPB 

1?  carbon  black 

A-172 

51.6? 

15u 

22.2?  2'J0p 

25?  CTPB 

1?  carbon  black 

that  the  meaa  and  extinguishment  behavior  of  the  PU-AP  propellant  is 
caused  by  an  accumulation  of  molten  binder  on  the  surface. 

A  number  of  analytical  modela  have  been  proposed  for  steady-state 
burning  of  solid  propellanta  (Ref.  24—2 7 »  see  also  references  cited  in 
Ref.  25).  Of  theae  the  beat  correlation  with  experimental  data  la  pro¬ 
vided  by  the  granular  diffusion  flame  (CDF)  model  proposed  originally  by 
Sunmerfield,  et  al  (Ref.  27). 
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FIG.  3.2.  Burning  Rate  Curve  for  a  Polyurethane 
Ammonium  Perchlorate  Propellant  Containing  1Z 
n-butyl  Ferrocene  Catalyst. 


The  GDF  model  predicts  that  the  burning  rate  will  be  related  to 
pressure  according  to  the  expression  p/r  *»  e  +  bp2'  .  Thia  expression 
predicts  that  plots  of  p/r  versus  p2' 3  will  be  straight  linea,  with  a 
single  value  of  the  intercept  &  for  a  given  chemical  composition  of  pro¬ 
pellant.  The  data  of  this  study,  when  plotted  in  this  manner,  give  the 
plota  of  Fig.  3.4  and  3.5. 


pressure  according  to  the  expression  p/r 
predicts  that  plots  of  p/r  versus  p2'3  wi 


As  can  be  seen,  the  GDF  model  does  not  implicitly  correlate  either 
set  of  data,  although  it  provides  a  rough  fit  to  the  CTPB  data.  The  dis¬ 
crepancy  between  the  CDF  model  and  the  PU-AP  data  is  drastic  and  is 
presumably  due  in  part  to  the  role  of  the  molten  binder.  The  discrepancy 
in  the  case  of  the  CTPB-AP  data  has  qualitative  importance  and  way  be 
due  to  the  blmodal  oxidizer  particle  size  distribution.  Particle  size 
distribution,  an  attribute  of  the  propellant  not  embodied  in  existing 
models,  is  important  when  trying  to  predict  the  burning  ratea  of  in- 
service  propellants  which  may  have  blmodal  or  triraodel  AP  particle  size 
distributiona.  In  b^th  cases,  the  microscopic  heterogeneity  of  the 
propellant  appears  to  have  a  substantial  effect  on  the  burning  rate 
characteristics,  but  is  absent  from  published  models. 

In  view  of  the  complicated  dependence  of  the  burning  rata  on  pres¬ 
sure  observed  in  this  study,  the  use  of  r  ■  cpn  or  p/r  »  a  +  bp2'3  for 
anything  other  than  a  first  approximation  is  unwise.  In  particular, 
the  use  of  r  -  cpn  in  the  several  perturbation  models  for  solid  propel¬ 
lant  combustion  inetablllty  must  lead  to  some  misrepresentation  of  the 
pressure  dependence  of  the  response  function. 

Further  details  on  this  work  are  available  in  Ref.  6. 


FIC.  3.3.  Burning  Rate  Curve*  for  Carboxyl-Te 
Polybutadlene-Aamonlum  Perchlorate  Propellant* 
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FIG.  3. A.  Data  for  the  Polyurethane-Ammonium  Perchlorate  Pro¬ 
pellants  Plotted  In  the  Manner  Suggested  by  the  Granular  Diffu¬ 
sion  Flame  Model  Correlation. 
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FIG.  3.5.  Data  for  the  Carboxyl-Terminated  Polybutadiene” 
Ammonium  Perchlorate  Propellanta  .  lotted  in  the  Manner  Sug¬ 
gested  by  t’-'?  Granular  Diffuaion  Flame  Model  Correlation. 
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3.3.  SURFACE  STRUCTURE  OF  COMPOSITE  SOLID  PROPELLANTS 

In  order  to  quantitatively  predict  the  combustion  behavior  of  solid 
propellants  it  is  necessary  to  reduce  the  description  of  this  behavior 
to  a  set  of  relations  which  are  amenable  to  mathematical  expression. 

The  surface  structure  of  a  burning  solid  propellant  must  be  considered 
when  attempting  to  mathematically  model  the  propellant  combustion.  The 
majority  of  models  proposed  for  ateady  and  nonateady  combustion  of  solid 
.propellants  treat  the  burning  aurface  as  s  planar,  dry,  homogeneous  en¬ 
tity  with  one  simple  Arrhenius  expressiou  describing  the  kinetica  of  the 
procesaea  over  the  entire  aurface.  For  composite  propellants,  thia  ob¬ 
viously  is  not  the  case.  One  justification  for  auch  simplification  ia 
that  the  complex  reality  ia  approximated  by  a  more  tractable  mathemati¬ 
cal  representation.  Another  reason  ia  that  detailed  Information  defin¬ 
ing  the  structure  of  the  surfsce  accurately  enough  to  construct  such  a 
model  is  not  available.  Because  of  the  latter  reaaon,  we  initiated  an 
experimental  investigation  to  provide  a  better  understanding  of  the 
physical  nature  of  the  aurface  and  determine,  if  possible,  the  extent 
of  heterogeneous  or  aubaurface  reactions.  In  this  atudy  burning  pro¬ 
pellant  samples  were  quenched  and  then  examined  using  an  SEM  to  infer 
the  structure  of  the  burning  aurfsce  prior  to  quench. 

The  unmetallized  propellants  rested  consisted  of  a  crystalline 
oxidizer  and  either  polyurethene  or  csrboxy-terminated  polybutadiene 
(CTPB)  binder.  The  propellants  were  formulated  with  either  s  unimodal 
oxidizer  particle  size  distribution  or  a  narrow-cut  bimodal  distribution, 
with  no  additives  such  as  burning  rate  cat  lysta  or  powdered  metal. 

Small  atranda  of  these  propellants  were  burned  si  pressures  between  one 
atmosphere  and  80C  psia  in  a  combustion  bomb  containing  an  inert  gaa. 

Once  ateady-atate  conditions  were  attained,  the  sample  waa  extinguished 
by  rapid  depresaurlzation  of  the  combustion  bomb.  The  extinguished 
samples  were  cut  to  an  appropriate  aize,  plated  with  a  gold-palladium 
alloy  or  a  few  ar^stroms  thickness,  and  examined  and  photographed  by 
using  the  SEM. 

Results  typical  of  this  study  are  shown  in  Fig.  3.6-3.10.  Figme 

3.6  presents  SEM  micrographs  of  aamplea  consisting  of  AP  and  a  polyure- 
tnane  binder  which  were  burned  at  pressures  of  100  psia  (Fig.  3.6a-b) 
and  200  psia  (Fig.  3.6c-f).  Inspection  of  the  photographs  reveala  that 
at  low  preaaurea  the  pyrolysis  rate  of  AP  waa  less  than  that  of  the 
polyurethane  binder,  thus  confirming  the  observation  reported  by  Baatress 
(Ref.  24)  that  the  AP  protrudea  above  the  binder  at  low  pressures.  Alao, 
the  smooth  appearance  of  the  binder  suggests  that  it  was  molten  during 
burning.  It  is  of  interest  to  note  thst  sll  of  the  views  shown  in 

Fig.  3.6  (PU-AP  propellant)  Indicate  definite  undercutting  of  the  binder- 
oxidizer  particle  boundary.  Such  undercutting  illustrates  the  error  of 
using  one  space  coordinate  to  try  to  describe  the  combustion.  Figure 

3.7  illustrates  the  effects  of  pressure  on  the  surface  structure.  The 
photogrspha  on  the  left  (Fig.  3.7a,  c  and  e)  are  of  samples  obtained  at 
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FTG.  3.6.  Scanning  Electron  Microscope  Micrographs 
of  Polyurethane-Ammonium  Perchlorate  Propellant 
Which  was  Quenched  While  Burning  at  IDO  pala  (a-b) 
and  200  pala  (c-f) . 


100  pala  while  those  S' "plea  on  the  right  (Fig.  3.7b,  d  and  f)  were  ob¬ 
tained  at  800  pala.  Both  aeta  of  samples  conclusively  demonstrate  that 
the  binder  was  molten  prior  to  quench.  At  the  lower  pressure  the  AP 
crystals  protrude  above  the  binder  aurfacc  and  undercutting  la  again 
noted.  At  the  higher  preaaure,  the  crystals  are  recessed  below  the 
binder  surface  to  the  extent  that  the  molten  binder  has  apparently  flowed 
into  the  recesses  and  covered  the  crystal  surfaces  (possibly  during  ex¬ 
tinguishment).  Thia  phenomenon  provides  cn  explanation  of  why  some 
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FIG.  3.7.  Scanning  Electron  Microscope  Micrographa 
of  Polyurethane-Ammonium  Perchlorate  Propellant 
'hich  was  Quenched  While  Burning  at  100  psia  (a,  c 
and  e)  and  800  psia  (b,  d  and  f )  . 

polyurethane  propellants  do  not  sustain  combustion  at  pressures  above 
approximately  600  psia  as  previously  reported  (Ref.  28).  As  the  pres¬ 
sure  is  increased  the  regression  rate  of  the  AP  eventually  becomes 
greater  than  the  binder  regression  and  the  molten  binder  then  flows  over 
the  AP  crystals  inhibiting  their  combustion  and  causing  extinguishment. 
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Earlier  we  showed  that,  under  certain  conditions,  AP  crystals  have 
a  surface  melt.  Prior  to  the  SOI  studies  described  here,  there  had 
been  considerable  controversy  surrounding  the  question  of  whether  these 
crystals  in  the  propellant  matrix  melted  or  sublimed  during  combustion. 
The  SEM  permitted  us  to  make  a  detailed  examination  of  quenched  propel¬ 
lant  samples  which  showed  that  the  AP  in  a  burning  propellsnt  wss  molten 
(Fig.  3.8  and  3.9).  The  bubble  formation  shown  in  Fig.  3.8  is  typical 
of  the  structure  seen  on  quenched  samples  of  AP  which  was  interpreted 


FIG.  3.8.  Scanning  Electron  Microscope  Micrographs 
of  Polyurethane-Ammonium  Perchlorate  Propellant 
Which  was  Quenched  While  Burning  at  100  psia. 


Iss  evidence  that  the  surface  of  the  AP  consists  of  a  thin  molten  layer 
during  Durning.  Gases  formed  from  the  decomposition  of  the  crystals 
within  the  molten  layer  expanded  during  rapid  depressurization  causing 
the  bubble  formation  as  the  melt  froze.  Indeed  it  wculd  be  difficult 
to  explain  the  bubble-like  structure  without  assuming  the  existence  of 
s  liquid  state.  Similar  types  of  structures  sre  evident  in  Fig.  3.6 
and  3.9  as  well  ss  in  the  majority  of  the  samples  studied.  The  volcar~- 
like  structure  seen  in  Fig.  3.6a,  e  and  f  and  Fig.  3.9,  and  the  vented 
structure  shown  in  Fig.  3.6c-d  also  indlcste  thst  subsurface  reactions 
have  taken  place  within  the  crystal  prior  to  quench.  These  structures 
are  in  agreement  with  those  observed  when  pure  single  crystals  of  AP  are 
burned,  quenched,  snd  examined. 
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FIG.  3.9.  Scanning  Electron  Microscope  Micrographs 
of  Carboxyl  Terminated  Polybutadiene-Ammonium  Per¬ 
chlorate  Propellant  Which  was  Quenched  While  Burning 
at  600  psia. 
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FIG.  3.10  Scanning  Electron  Microscope  Microgrsph 
of  sn  HMX-Poly urethane  Propellant . 


By  mechanically  stressing  the  quenched  samples  it  was  possible  to 
break  the  AP-binder  bond,  releasing  the  crystals  as  seen  In  Fig.  3.8. 
Examination  of  the  displaced  crystals  shown  in  Fig.  3.8c  and  one  of  the 
crsters  left  undisturbed  in  the  binder  ss  seen  in  Fig.  3.8d  shows  no 
indication  of  subsurface  or  interfacial  reaction  between  the  AP  and  the 
binder. 

In  order  to  explore  the  effect  of  binder-type  on  the  surface  struc¬ 
ture,  s  propellant  containing  s  carboxy-tenninsted  polybutadiene  (CTPB) 
binder  wss  tested.  Results  from  these  tests  st  600  psls  are  included  in 
Fig.  3.9.  The  micrographs  indicate  that  the  surface  of  the  AP  crystals 
wss  molten  snd  thst  subsurface  reactions  in  the  molten  phsse  of  the  AP 
resulted  in  gas  libers tion  within  the  molten  phsse.  The  microgrsph  seen 
in  Fig.  3.9d  snd  f  seems  to  indicate  thst  lnterfscisl  reactions  were  not 
present  st  the  binder-oxidizer  (CTPB-AP)  interface  st  this  pressure.  No 
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undermining  was  observed  on  samples  obtained  at  other  pressures,  indi¬ 
cating  a  lack  of  interfacial  reactions  with  the  CTPB  samples.  Also, 
there  is  no  indication  of  the  CTPB  binder  melting  as  with  the  polyure¬ 
thane  binder  (i.e.,  the  binder  does  not  have  a  smooth  appearance  nor  the 
"cracked-mud"  appearance  of  Fig.  3.7). 

Study  of  the  behavior  of  oxidizers  other  than  AP  in  the  combustion 
zone  has  been  rather  limited.  However  the  various  crystalline  oxidizers 
of  intereat  differ  widely  in  their  high  temperature  behavior,  and  an  ob¬ 
jective  view  of  combustion  zone  processes  cannot  be  formulated  on  the 
basia  of  AP  propellants  alone.  Figure  3.10  shows  SEM  pictures  of  the 
quenched  surface  of  a  propellant  with  a  polyurethane  binder  and  HMX 
(cyclotetramethylenetetranitramine) .  The  melting  point  of  HMX  ia  ao 
low  (280°C)  that  it  would  be  molten  to  an  appreciable  depth  beneath  the 
propellant  surface.  From  examination  of  Fig.  3.10  it  is  evident  that 
tbe  KMX  was  decomposing  at  an  appreciable  depth  below  the  aurface,  with 
a  very  "turbulent"  surface  (Fig.  3.10a-b)  and  bubbles  to  considerable 
depth  (in  cross-sections  of  the  sample.  Fig.  3.10c-d). 

These  investigations  of  surface  structure  of  the  combustion  zone 
are  continuing,  including  cooperation  with  an  Air  Force  program  at  the 
Lockheed  Propulsion  Co.  Further  details  of  the  work  may  be  obtained 
from  Ref.  16,  and  29-31. 


3.4.  METAL  AGGLOMERATION 

The  utility  of  including  metallic  fuels  in  propellants  depends 
upon  how  efficiently  the  metal  ignites,  burns,  and  produces  high  tem¬ 
perature  working  gases.  Aluminum  as  a  fuel  is  being  studied  in  detail 
because  it  is  in  widespread  use,  it  frequently  exhibits  poor  combustion 
efficiency  and  it  is  an  important  factor  in  driving  axial  mode  and  low 
frequency  combustion  instability.  The  current  work  is  an  investigation 
of  the  preignition  behavior  of  aluminum,  using  a  hot  stage  microscope, 
a  gas  laser,  and  the  agglomeration,  ignition,  and  combustion  of  aluminum 
on  and  near  the  aurface  of  composite  propellants. 

3.4.1.  Particle  Prelgnitlon  Studies 

The  experimental  setup  for  the  preignition  studies  of  aluminum 
particles  is  described  in  Ref.  32.  Two  heating  methods  have  been 
applied,  namely  heating  on  a  microscope  hot  stage  and  heating  by  laser 
radiation.  The  first  method  produces  only  low  heating  rates  (<900°C/min) 
but  the  tempersture  can  be  readily  measured  whereas  the  laser  heating 
yields  higher  heating  rates  but  makes  temperature  measurements  very 
difficult  or  even  Impossible.  After  heating  and  cooling,  the  samples 
were  examined  with  the  optical  microscope  as  well  as  the  SEM.  The  re¬ 
sults  are  presented  in  detail  in  Ref,  34  and  will  be  described  here  only 
briefly. 


25 


NWC  TP  4749 


Heating  of  spherical  aluminum  particles  (105-125p  diameter)  in  air 
at  the  relatively  low  rate  of  about  900°C/min  produced  surface  oxidation; 
in  the  areas  where  per  tides  were  touching,  agglomeration  at  temperatures 
as  low  as  the  melting  point  of  aluminum  (Tm  *  659. 7°C)  occurred.  At  a 
constant  heating  rate  the  amount  of  oxidation  is  a  function  of  temper¬ 
ature,  as  shorn  in  Fig.  3.11.  Figure  3.11a  shows  the  smooth  but  wrinkled 
surface  of  two  agglomerated  particlea  after  heating  to  690°C.  Heating  to 
845°C  did  not  change  the  smoothness  of  the  oxide  akin  but  resulted  in 
severe  cracking  and  wrinkling.  This  can  be  aeen  in  Fig.  3.11b  of  a 
"bridge"  connecting  two  agglomerated  particles .  In  both  caaes  wrinkling 
waa  observed  to  occur  during  cooling  of  the  particles  and  hence  is  not 
of  great  concern  to  the  preignition  behavior.  At  1040°C  the  particles 
were  oxidized  to  the  extent  that  the  resulting  oxide  sheila  were  rigid 
and  no  longer  could  wrinkle  on  cooling.  Consequently,  the  particlea  re¬ 
tained  their  spherical  shape  and  exhibited  a  coarse  granular  structure 
of  oxide  surface  (Fig.  3.11c).  Some  large  spherical  agglomerates  were 
obtained.  The  spherical  shape  indicated  that  solid  oxide  (melting  point 
Tm  -  2042°C)  was  incorporated  in  the  agglomeratea.  Finally,  heating  to 
1200°C  increased  oxidation  to  the  point  where  the  oxide  on  the  agglomer¬ 
ate  surface  was  granular  (Fig.  3. lid)  and  much  coaraer  than  that  ob¬ 
tained  at  1040°C. 

Laser  heating  of  spherical  aluminum  particles  in  air  at  high  rates 
to  temperatures  below  the  melting  point  of  aluminum  oxide  (T  <  2042°C) 
resulted  in  heavy  oxidation  and  agglomeration  (Fig.  3.12a)  very  similar 
to  that  obtained  on  the  hot  stage  microscope.  Laser  heating  to  temper¬ 
atures  above  the  melting  point  of  aluminum  oxide  (T  >  2042°C)  led  to 
fa3t  oxidation  of  the  individual  metal  particles  with  greatly  reduced 
agglomeration  (Fig.  3.12b).  Heavy  smoke  indicated  the  transition  to 
burning. 

In  summary,  agglomeration  on  the  hot  stage  occurred  at  rather  low 
temperatures  when  oxidation  was  atlll  at  low  level  and  the  particles 
were  in  physical  contact.  Below  845°C  little  oxidation  occurred  which 
favored  agglomeration  in  that  temperature  range.  Heavy  oxidation  made 
the  oxide  ahella  rigid  and  reduced  agglomeration  at  high  temperatures. 
Heating  rates  (by  laaer  radiation)  above  the  melting  point  of  aluminum 
oxide  favored  the  oxidation  of  the  individual  aluminum  particle.  No 
large  agglomerates  involving  the  flow  of  melted  aluminum  were  observed. 
Instead,  the  oxide  seemed  to  fuse  together  forming  an  Irregularly  shaped 
structure. 

3.4.2.  Aluminum  Behavior  in  Propellants 

In  preceding  reports  a  qualitative  description  of  aluminum  accumu¬ 
lation  on  the  surface  of  burning  propellants  (referred  to  as  the  "pocket 
model")  (Ref.  4)  and  the  qualitative  description  of  the  microscopic  de¬ 
tails  of  the  agglomeration  processea  (Ref.  5)  were  presented.  Subse¬ 
quent  work  in  this  area  haa  been  intermittent  due  to  reduced  funding; 
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FIG.  3.12.  Aluminum  Particles  After  Heating  in  Air 
hy  a  Laser  Beam  to  Temperatures,  (a)  Below  and  (b) 
niove  Melting  Point  of  Aluminum  Oxide. 


it  has  been  aimed  at  (1)  obtaining  more  conclusive  microscopic  evidence 
of  the  sluminum  sintering  process,  (2)  measurement  of  aluminum  agglomer¬ 
ate  sizes  from  high  speed,  magnified  motion  pictures  to  provide  a  more 
quantitative  basis  for  the  "pocket  model",  and  (3)  combining  all  past 
work  into  a  comprehensive  open  literature  publication. 

Earlier  attempts  to  obtain  partially  burned  samples  of  propellsnt 
for  microscopic  examination  of  the  aluminum  agglomerates  (Ref.  5  and  33) 
by  rapid  depressurization  and  by  thermal  quenching  on  a  stainless  steel 
plate  resulted  in  samples  which  did  not  show  the  details  of  the  sinter¬ 
ing  process  ss  conclusively  sa  desired.  Subsequently,  samples  of  alumi¬ 
nized  propellants  were  burned  and  thermally  quenched  using  a  copper  vise 
heat  sink  (Ref.  7).  Scanning  electron  microscope  pictures  of  the  sam¬ 
ples  quenched  in  this  way  are  shown  in  Fig.  3.13  and  3.14;  for  comparison. 
Fig.  3.15  shows  SEM  photographs  of  100  micron  sluminum  balls  sintered 
together  by  Impingement  on  s  sapphire  piste  hested  to  1400°C  (Ref.  5). 


Figure  3.13  shows  successively  magnified  pictures  from  the  SEM  of  one 
ares  of  the  burning  surface  of  thermally  quenched  A-49  propellsnt  which  con¬ 
tained  16X  15  micron  sluminum,  59X  200  micron  AP,  and  25X  polyurethane  binder. 
An  AF  particle  approximately  one-hslf  consumed  is  evident  in  the  upper  center 
of  Fig.  3.13s;  the  accumulation  on  top  of  the  partially  consumed  AF  is 
aluminum  particles  of  approximately  15  micron  diameter  thst  sre  probably 
sintered  together.  The  accumulate  ss  s  whole  appears  to  be  held  onto  the 
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FIG.  3.13.  Surface  of  Thermally  Quenched  Propel¬ 
lant  Shoving  Accumulation  of  Aluminum  of  the  Surface 
of  the  AP  Particles  and  Binder. 
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aurface  of  the  AP  particle  by  being  vetted  by  the  liquid  AP  surface 
Ref.  13).  In  Fig.  3.13b  alntered  aluminum  particles  can  be  aeen  on  the 
binder  surface  (in  the  foreground  of  the  AP  particle) ;  these  particlea 
also  appear  to  be  held  on  the  binder  aurface  by  being  wetted  by  a  liquid 
binder  surface.  Thus,  these  pictures  are  rather  conclusive  evidence 
that  a  major  portion  of  the  aluminum  la  retained  on  the  aurface  of  the 
propellant,  either  on  the  "wet"  binder  or  'Vet"  AP  aurface.  It  appeara 
that  a  very  important  factor  in  the  agglomeration  proceaa  is  that  dur¬ 
ing  the  time  between  the  emergence  of  an  aluminum  particle  on  the  aurface 
and  lta  inclusion  as  part  of  an  agglomerate.  It  is  held  on  the  aurface 
by  thia  "wetted"  condition.  As  the  surface  continues  to  recede,  the  pop¬ 
ulation  of  the  accumulated  aluminum  Increases  until  the  particles  are  in 
contact  with  one  another.  During  this  accumulation  time,  oxidative  aelf- 
heating  and  increased  proximity  to  the  gaa  phase  heat  aource  ralae  the 
temperature  of  the  aluminum  to  its  melting  point;  at  the  melting  point 
of  aluminum,  the  sintering  or  "bridging"  of  the  aluminum  particlea  begins 
(Ref.  5). 
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Figure  3.14  shows  SEM  photographs  of  s  thermally  quenched  sample  of 
A— 107  propellant  which  contained  16%  15  micron  aluminum  snd  59%  600  micron 
AP.  An  sssemhlsge  of  loosely  sintered  aluminum  particles  is  shown  in 
Fig.  3.14a.  Successive  magnifications  sre  focused  on  the  "bridge"  of 
aluminum  snd/or  aluminum  oxide  connecting  two  aluminum  particles.  This 
series  of  pictures  when  compared  to  those  in  Fig.  3.15  show  the  geometri¬ 
es  1  similarity  between  the  "bridging"  process  on  the  surface  of  burning 
propellants  snd  the  more  easily  Interpreted  "bridging"  process  occurring 
during  the  rapid  heating  of  100  micron  aluminum  spheres  in  sir  (Ref.  5). 
Thus,  the  pictures  in  Fig.  3.24  sre  rather  conclusive  evidence  that  the 
sintering  process  thst  leads  to  formstion  of  aluminum  agglomerates  in 
burning  propellants  is  very  closely  related  to  the  mechanistic  descrip¬ 
tion  of  the  aluminum  "bridging”  process  described  in  Ref.  5. 


FIG.  3.14.  Aluminum  Agglomeration  on  Surface  of 
Propellant  Showing  the  "Bridging”  Together  of 
Particles. 


o  b 

FIG.  3.15.  Aluminum  "Bridging  of 
Heated  Single  Particles  (Ref.  5). 


During  the  course  of  organizing  the  material  on  aluminum  agglomer¬ 
ation  for  an  open  literature  publication,  it  became  obvious  that  de¬ 
tailed  measurement  of  the  agglomerate  size  distributions  was  very  de¬ 
sirable.  Therefore,  motion  picture  sequences  from  18  tests  of  14 
different  propellant  formulations  were  selected  for  measurement.  Most 
of  the  films  have  been  measured  and  although  there  appears  to  be  aome 
discrepancies  in  the  data,  it  appears  that  the  quantitative  information 
from  the  measurements  ia,  in  general,  conaiatent  with  the  qualitative 
picture  of  the  "pocket  model". 

The  procedure  for  determining  the  size-number  distribution  from 
the  movie  film  was  to  count  each  agglomerate  and  measure  lta  apparent 
diameter  in  a  selected  volume  at  some  arbitrary  time  and  to  continue  to 
count  and  measure  the  additional  incoming  agglomerates  in  subsequent 
consecutive  frames  until  a  minimum  of  200  agglomerates  had  b.*en  measured. 
The  volume  selected  for  measurement  was  4000  microns  wide  (slightly 
narrower  than  the  film  image  width),  \1500  microns  thick  (the  thickness 
of  the  original  propellant  sample) ,  and  1000  microns  above  the  burning 
surface.  The  apparent  agglomerate  diameter  was  recorded  in  10  micron 
increments.  Representative  data  are  presented  in  Fig.  3.16  where  the 
ordinate, N' ,  of  the  histogram  la  the  number  of  agglomerates  normalized 
by  the  count  of  the  moat  numerous  size. 

Since  the  "pocket  model"  states  that  the  size  of  aluminum  agglomer¬ 
ate  is  principally  determined  by  the  total  volume  of  aluminum  contained 
in  the  binder  pocket  between  the  .VP  particles  and  since  the  behavior  of 
the  totsl  mass  of  aluminum  ia  the  quantity  of  intereat,  it  is  necessary 
to  convert  the  number-diameter  relationship  to  a  relationship  between 
the  maas  of  the  agglomerate  and  the  size  distribution.  The  data  in 
Fig.  3.16  were  converted  to  reflect  the  role  of  the  mass  of  the  agglomer¬ 
ates  as  shown  in  Fig.  3.17.  In  this  figure,  M’  is  the  normalized  total 
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mass  of  the  aluminum  agglomerate;  it  was  obtained  by  multiplying  the 
normalized  number  of  agglomerates,  N',  by  the  cuba  of  the  appropriate 
diameter  and  normalizing  with  respect  to  the  largest  mass.  In  order  to 
examine  what  effect  a  particular  propellsnt  vavisble  had  on  aluminum 
agglomeration,  it  was  necessary  to  represent  the  agglomerate  mass  dis¬ 
tribution  by  a  single  number.  The  agglomerate  diameter  corresponding 
to  the  average  normalized  mass  of  all  the  agglomerates,  i.e.,  the 
centroid  with  respect  to  the  diameter  of  the  M'  versus  agglomerate  diam¬ 
eter  histogram,  was  chcsen  as  the  variable  for  describing  the  agglomerate 
mass.  Work  Is  in  progress  on  reducing  the  data  to  the  form  of  "mass 
sverage  diameter”  as  a  function  of  propellant  parameters. 

An  example  of  the  type  of  information  being  obtained  is  shown  in 
Fig.  3.18  where  the  "mass  average  diameter",  d',  is  plotted  against  the 
AP  particle  size.  The  trend  of  increasing  agglomerate  size  with  in¬ 
creasing  AP  size  is  in  agreement  with  the  pocket  model.  Similar  data 
is  being  obtained  as  a  function  of  original  aluminum  particle  size, 
percentage  of  aluminum,  pressure,  binder  percentage  and  bindci  type. 

A  paper  for  the  opeti  literature  covering  all  the  work  on  aluminum 
behavior  is  being  prepared;  its  completion  date  is  contingent  on  the 
completion  and  reconciliation  of  the  reduced  data. 

3.4.3.  Interpretation  of  Results 

The  results  obtained  by  the  various  investigations  provide  a  quali¬ 
tative  understanding  of  how  and  why  aluminum  accumulates  and  agglomer¬ 
ates  in  the  combustion  zone.  Since  aluminum  is  a  pure  element  with  a 
high  boiling  point,  it  does  not  decompose  and/or  gasify  at  the  propel¬ 
lant  surface  temperatures  (''•600°C)  .  It  comes  to  the  surface  with  a  high¬ 
ly  protective  coating  of  oxide,  and  is  surrounded  by  a  flow  of  fuel 
vapors  from  the  binder.  Conditions  are  not  particularly  favorable  for 
oxidation,  although  It  may  be  appreciable  when  the  aluminum  is  present 
in  very  fine  particle  size  (e.g.,  <  5  p) .  Since  the  aluminum  is  seen  to 
accumulate  on  the  surface,  there  is  apparently  enough  surface  force  to 
retain  emerging  particles  (one  would  expect  this  to  be  dependent  on  the 
type  of  binder) . 

As  the  aluminum  particles  sccumulate  on  the  surface,  they  rise 
further  in  temperature  due  to  protrusion  towards  the  diffusion  flame,  and 
they  crowd  together  as  more  particles  emerge  from  the  receding  binder. 

From  studies  described  above.  It  seems  clear  that  the  particles  will 
"sinter"  together  when  the  aluminum  melting  point  is  reached,  this  process 
taking  place  by  oxidation  of  aluminum  leaking  through  the  expansion 
cracks  in  the  oxide  shell.  This  leads  to  an  Irregular  array  of  loosely 
held  alumluum  on  the  surface,  which  may  continue  to  oxidize  as  the 
■proximity  to  the  diffusion  flame  Increases  and  the  oxide  continues  to 
crack  due  to  difference  in  expansion  coefficient  of  A1  and  AI2O3.  From 
motion  pictures  and  quenched  samples,  it  is  clear  that  the  metal  resides 
for  an  appreciable  time  on  the  surface  without  inflammation.  The  material 
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is  not  mobile  in  thia  phase  of  the  agglomeration.  The  increasing  obstruc¬ 
tion  of  gas  flow  by  the  growing  agglomerate  eventually  provides  forces 
sufficient  for  detachment  from  the  surface.  This  trend  ia  accompanied 
by  chemical  changes  in  the  underlying  material  sud  surrounding  gas  flow, 
due  to  consumption  of  binder  filaments  and  oxidizer  particles.  It  would 
be  difficult  to  know  to  what  extent  separation  from  the  aurface  la  due 
to  lncreaaing  obstruction  of  flow,  and  to  what  extent  due  to  changing 
surface  forces  and  chemical  environment.  It  aeems  likely  that  the  con¬ 
sumption  of  a  volume  element  of  binder  between  oxidizer  particles  will 
leave  an  accumulate  that  consists  of  most  of  the  metal  from  that  fuel 
pocket.  The  emergeuce  of  an  underlying  oxidizer  particle  provides  a 
transition  in  both  fluid  dynamic  and  chemical  environment  conducive  to 
inflammation.  This  hsa  lead  to  the  "pocket”  model  mechanism  for  metal 
agglomeration,  and  the  results  of  the  present  studies  on  propellants 
with  unlmodsl  oxidizer  particle  size  distribution  (which  led  to  the 
"pocket"  model)  indicate  that  the  foregoing  microstructure-related  argu¬ 
ment  for  agglomeration  is  relevant.  However,  results  with  propellants 
having  polymodal  oxidizer  particle  size  (and  no  real  "pockets")  show 
agglomeration  to  the  same  range  of  droplet  sizes,  so  that  other  factors 
such  as  protrusion  of  agglomerates  into  the  diffusion  flame  can  also  be 
importsnt  fsctora  in  the  agglomeration  sequence.  Further,  at  low  pres¬ 
sures  the  agglomeration  csn  be  much  more  extensive  because  of  the  leas 
favorable  conditions  for  inflammation,  and  agglomerate  droplets  may  re¬ 
sult  from  inflammation  of  fsirly  extensive  sintered  layers  of  aluminum. 
This  procesa  is  particularly  conspicuous  in  "propellants"  made  by  press¬ 
ing  powdered  AF  and  aluminum  without  binder.  Pellets  made  in  this  way 
exhibit  massive  agglomeration  of  aluminum,  with  the  sintered  layer  being 
consumed  intermittently  by  spread  cf  inflammation  from  local  initiation 
altea. 

In  most  situations  the  agglomerate  la  seen  to  ignite  in  the  process 
of  leaving  the  surface,  so  thst  inflammation  and  separation  constitute 
one  complicated  event.  The  aintered  accumulate  Is  seen  to  draw  (or  roll) 
up  into  a  burning  agglomerate  droplet,  which  usually  then  leaves  the 
surface  (window  bomb  tests),  but  in  some  situations  Instead  rolls  around 
picking  up  more  accumulate.  Under  some  conditions  particularly  unfavor¬ 
able  to  aluminum  combustion,  the  agglomerate  does  not  become  fully  in¬ 
flamed,  but  instead  remains  a  lumlnoua,  irregular  agglomerate  which  hsngs 
tenuously  to  the  surface  and  eventually  flosts  sway  without  becoming  a 
droplet. 

Inflammation  ia  sn  easily  distinguishable  event  in  most  situations, 
because  it  involves  the  transition  from  a  slowly  reacting,  solid  oxide- 
costed  agglomerate  to  a  rapidly  reacting  droplet  In  which  the  oxide  has 
melted  and  partly  retracted  from  much  of  the  aluminum  surface.  This 
stage  is  easily  visible  because  the  sgglomerste  becomes  aphericsl,  ex¬ 
tremely  luminous,  has  s  visible  oxide  amoke  trail,  and  the  surface  is 
clearly  comprised  of  aresa  of  oxide  (high  emissivity)  and  metal  (low 
emiaaivity,  with  reflected  highlights  under  external  illumination). 

Theae  observations  are  similar  to  those  made  in  independent  studies  of 
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combustion  of  single  alumioum  particles  io  controlled  stmospheric  ex- 
perimeots  (Ref.  32,  34  aod  35)  where  oxygeo  is  the  priocipsl  oxidizing 
specie. 

The  mobility  of  the  sgglomerste  incresses  duriog  or  after  inflam¬ 
mation:  in  the  flow  field  of  the  rocket  motor  the  droplet  probably 
lesves  the  surface  st  this  poiot,  but  io  window  bomb  tests  it  sometimes 
wsnders  over  the  surface  with  further  sgglomerstioo.  From  the  stsodpoiot 
of  oversll  combustion  behavior,  this  poiot  is  probably  quite  important. 
Further  sgglomerstion  cao  sffect  efficieocy  of  combustion  of  the  metal 
io  the  motor.  Retention  of  burning  metsl  oesr  the  surface  probsbly 
affects  the  burniog  rste  of  the  propellsot  (io  s  pressure-dependeot  way). 
Susceptibility  of  this  loosely  held  sgglomerste  to  premature  ioflammation 
snd/or  sepsration  from  the  surface  by  s  flow  disturbance  may  cause  un¬ 
stable  burning.  All  of  these  factors  sre  oo  doubt  highly  dependent  oo 
such  vsrisbles  as  metsl  particle  size,  oxidizer  particle  size,  stoichio¬ 
metry,  binder  melting  and  decomposition  characteristics,  sod  gas  flow 
environment.  Tecsuse  of  the  lsrge  number  of  vsrisbles  iovolved,  the 
continuing  ignorance  of  binder  sod  oxidizer  pyrolysis,  sod  the  limited 
knowledge  of  the  diffusioo  flsme  structure,  the  cootribution  of  the 
metsl  to  the  vsrious  sspects  of  propellsnt  combustloo  remains  obscure. 
However,  increasing  quslitstive  understanding  holds  promise  for  rational 
ordering  of  experiments!  observations  sod  for  development  of  techniques 
to  cootrol  metsl  combustion  behavior. 
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4.  BULK  MODE  INSTABILITY 


Oscillstory  combustion  in  low  frequency,  non-wave  modes  is  a 
problem  confined  primarily  to  motors  operating  at  low  combustor  L*, 
typical  of  spsce  motors  and  thrust  magnitude  control  motors  operating 
in  low  thrust  mode.  Analytical  and  laboratory  testing  of  this  phenome¬ 
non  has  been  relatively  successful  because  it  is  simpler  than  wave  mode 
Instability.  A  number  of  papers  have  been  published  recently  on  this 
problem,  snd  recent  work  on  this  program  is  summarized  in  some  of  these 
references  (Ref.  36-38),  and  will  not  be  repeated  here.  Instead,  a 
diagnosis  of  the  state  of  knowledge  and  future  needs  is  presented.  This 
diagnosis  emerged  from  a  critical  review  of  work  of  all  agencies  to  date, 
s  review  that  was  stimulated  by  a  judgement  that  the  experimental  results 
were  telling  the  scientific  community  much  more  thsn  we  were  hearing, 
and  that  our  thinking  was  unduly  inhibited  by  reliance  on  overly  simpli¬ 
fied  analytical  models. 

4.1.  STATUS  OF  KNOWLEDGE 

Analysis  of  BMI  consists  of  an  application  of  conservation  laws  to 
the  mass  and  energy  transport  in  the  combustor.  Modeling  of  the  combus¬ 
tion  zcne  itself  (usually  assumed  to  be  localized  at  the  propellant  sur¬ 
face)  is  the  object  of  an  additional  collection  of  studies  (Ref.  4  and  39) 
and  may  be  incorporated  into  the  combustor  analysis  by  use  of  the 
"response  function"  concept,  which  represents  the  combustion  zone  in 
terms  of  its  acoustic  response,  i.e.,  amplitude  and  phase  of  the  burning 
rate  oscillation  produced  by  a  given  pressure  oscillation.  I, inear  analy¬ 
ses  of  mass  conservation  in  the  combustor  lead  to  various  forms  (Ref.  36 
and  40-45)  of  the  equations 


R  cos  u>t  —  r  —  <t>  (aL*/c*) 


(4.1) 


R  sin  (or  ■  4  (wL*/c*) 


(4.2) 


where  a  and  u  sre  the  growth  rate  and  frequency  of  the  pressure  oscilla¬ 
tions;  L*  describes  the  burner  geometry  (free  volume/throat  area);  R  and 
ut  are  the  magnitude  and  phase  of  the  combustion  oscillations;  c*  is  the 
usual  effective  discharge  velocity  of  the  nozzle;  and  r  and  $  have  values 
near  one,  and  are  functions  of  the  ratio  of  specific  heats  (exact  values 
slightly  sensitive  to  ratio  of  combustor  flushing  time  to  period  of  os¬ 
cillation).  The  qualitative  validity  of  these  equations  is  apparently 
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not  suapect,  except  when  the  volume  of  the  combuatlon  zone  la  compnrable 
to  the  combustor  volume  (e.g.,  with  metallized  propellants  at  low  com¬ 
bustor  pressure) ,  Application  of  the  equations  Is  of  course  dependent 
on  knowledge  of  R  and  t  as  functions  of  u  and  other  relevant  variables. 

Perhapa  the  largest  source  of  uncertainty  in  BMI  analyses  is  the 
limited  knowledge  of  the  response  function  and  the  qualitative  sensi¬ 
tivity  of  the  solutions  of  Eq.  4.1  and  4.2  to  this  range  of  uncertainty 
of  R  and  t.  The  qualitative  dependence  of  R  and  t  on  frequency  Is  aeen 
In  both  theory  and  experiment  tc  be  of  the  form  3hown  in  Fig.  4. In. 

The  corresponding  solutions  to  Eq.  4.1  and  4.2  are  shown  In  Fig.  4.1b. 

The  calculations  using  Eq.  4.1  and  4.2  may  proceed  In  either  direction 
from  Fig.  4.1a  to  4.1b  or  from  4.1b  to  4.1a,  according  to  the  needs  of 
the  situation  (Ref.  36).  (In  the  example,  the  curves  in  Fig.  4.1b  were 
proposed  first  on  tie  basis  of  a  speculative  interpretation  of  certain 
experimental  measurements  of  a,  w  and  L*  mentioned  later,  and  the  form 
of  the  corresponding  calculated  response  functions  was  used  as  a  test 
of  the  plausibility  of  the  postulated  a-L*  curves.)  The  reaults  suggest 
that  the  form  of  the  aolution  a  versus  L*  differs  qualitatively  for 
rather  modest  changea  in  the  response  function. 

Paat  studies  of  BMI  have  represented  the  response  function  in  terms 
of  one  of  the  "QSHOD"  analytical  models  for  combustion  perturbations 
(quasi-steady,  gas  phase,  homogeneous  propellant,  one-dimensional  models) 
(Ref.  39  and  46).  These  models  have  been  used  in  papers  on  BMI  pri¬ 
marily  to  analyze  the  stability  trends  defined  by  the  a  “  0  condition. 

But  they  can  be  used  also  to  develop  a  complete  solution  such  as  the 
normalized  a  versus  L*  functions  illustrated  in  Fig.  4.2.  These  func¬ 
tions  were  calculated  from  Eq.  4.1  and  4.2,  the  definitions  of  A  and  L* 
indicated  in  the  figure,  and  Eq.  17  and  43  of  Ref.  39  (the  quantities 
A  and  B  are  regarded  in  the  QSHOD  models  as  properties  of  the  propellant, 
relatively  insensitive  to  pressure  or  frequency,  see  Ref.  39).  From 
Fig.  4.2  it  is  evident  that  considerable  variation  in  the  form  of  the 
a  -  l*  function  may  be  expected  even  within  the  constraints  of  the  QSHOD 
models ,  a  point  that  has  apparently  escaped  previous  notice. 

While  the  adequacy  of  the  QSHOD  models  will  be  questioned  later 
(see  also  Ref.  11,  38,  46  and  47),  it  is  appropriate  to  consider  first 
the  completeness  and  validity  of  past  applications  using  those  models. 
Much  of  the  earlier  work  on  BMI  relates  to  the  effect  of  combustor  pres¬ 
sure.  Recalling  that  the  constants  A  and  B  in  the  models  are  relatively 
pressure  independent  (Ref.  39),  the  effect  of  pressure  enters  in  only 
through  the  mean  burning  rate,  Y,  In  the  variables  A,  L*  and  SI  (i.e., 
the  solutions  in  Fig.  4.2  are  each  valid  for  all  pressures  with  the 
different  curves  presumably  representing  different  propellanta) .  The 
detailed  behavior  as  a  function  of  pressure  can  then  be  examined  in  any 
desired  context  by  use  of  the  appropriate  curve  in  Fig.  4.2  and  the 
relevant  steady-state  burning  rate-pressure  dependence  of  the  propellant 
wherever  burning  rate  appears  in  A,  L*  and  W.  As  an  example,  the 
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FIG.  4.2.  Solutions  to  Eq.  4.1  snd  4.2  for  Response  Functions 
From  QSHOD  Combustion  Model  With  Various  Vslues  of  the  Parameter  A 
(Which  is  Related  to  the  Kinetics  of  the  Surface  Decomposition 
Reactions).  The  vslue  of  the  parameter  B  in  Eq.  5.1  used  here 
wss  0.8. 


stability  limit  a  «  0  corresponds  to  specific  vslues  of  0 
L*  -  K$L*/r2c*  for  sny  particular  curve  in  Fig.  4.2,  i.e. 


(c*L*  n/<t)r  «  (C2c*L*  n/ic«)p 


where  r  wss  replaced  on  the  right-hsnd  side  by  a  frequently  used  pres¬ 
sure  dependence  r  ■  Cpn  (see  comment,  Section  3.2).  It  is  these  two 
relations  that  have  been  used  to  correlate  experimentally  determined 
trends  of  BMI.  Thus  in_Ref.  40,  Eq.  4.4  wss  used  to  correlate  stability 
limit  dsts  in  terms  of  p  snd  L*,  while  dsts  relating  w  to  r  snd  p  were 
correlated  with  the  stsbility  limit  defined  by  Eq.  4.3  in  Ref.  36,  38  snd 
42.  This  is  s  rsther  superficial  treatment  of  both  theory  snd  experi¬ 
ment  ,  which  contain  much  more  of  vslue  than  the  functional  form  of  the 
a  »  0  stability  limit. 
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4.2.  TROUBLES  WITH  PRESENT  STATUS  OF  KNOWLEDGE 

There  are  several  aspects  of  the  present  situation  that  seem  to  be 
blocking  further  progress.  One  Is  the  preoccupstion  with  stability 
Halts.  Most  experiments  yield  measurements  of  a,  u,  L* ,  p  and  r ,  snd 
these  measurements  csn  be  correlated  in  a  variety  of  ways  thst  sre  more 
meaningful  than  mere  verification  of  a  stability  limit.  As  sn  example, 
values  of  R  snd  t  csn  be  calculated,  and  In  principle  s  portion  of  tbe 
response  function  curve  csn  be  determined.  Only  limited  efforts  slong 
these  lines  have  been  reported  (Ref.  36  snd  39).  Even  In  the  discus¬ 
sions  of  stsblllty  limits,  undue  rellsnce  Is  made  on  the  results  pro¬ 
vided  by  the  QSHOD  models  (Eq.  4.3  snd  4.4).  The  resulting  comparisons 
between  experiment  snd  analysis  sre  sometimes  quite  unsatisfactory 
(Ref.  11,  38  snd  46).  A  particularly  good  example  Is  provided  by  BMX 
studies  (Ref.  38)  of  two  simllsr  propellant  formulations  in  which  the 
principal  difference  was  between  the  polyurethane  and  CTPB  binders.  The 
stability  limit  trends  predicted  by  the  QSHOD  theory  were  in  agreement 
with  the  results  from  the  CTPB  propellant,  but  in  sharp  disagreement 
with  the  results  from  the  polyurethsne  propellant. 

A  second  problem  with  recent  literature  concerns  the  use  of  the 
criterion  o  »  0  to  define  the  stability  limit.  This  criterion  defines 
a  minimum  frequency  ft  below  which  a  <  0  so  that  oscillations  do  not 
grow.  It  can  be  seen  from  Fig.  4.1  snd  4.2  thst  the  maximum  vslue  of  L* 
for  oscillatory  behsvlor  sometimes  corresponds  to  3a/3L*  -  »  rather  than 
a  »  0.  This  point  spparently  has  not  been  noticed  in  previous  dis¬ 
cussions  of  stsblllty  limit,  although  the  behavior  when  L*  is  near  this 
limit  Is  quail tstively  different  from  the  o  -  0  esse  snd  may  explain  some 
of  the  dissgreement  between  theory  and  experiment  noted  in  earlier  work 
(Ref.  11,  38.  46  snd  47).  Within  the  QSHOD  "theory"  the  principal 
difference  is  the  finite  vslue  of  a  st  the  stability  limit  (the  pressure 
dependence  remains  as  in  Eq.  4.3  snd  4.4).  Outside  the  QSHOD  theory,  it 
may  happen  thst  s  gradual  transition  from  Type  s  to  Type  b  to  Type  c 
behavior  (Fig.  4.1b)  occurs  as  a  function  of  s  test  variable  such  ss 
pressure,  with  a  corresponding  shift  from  sn  a  ■  0  type  of  stsblllty 
limit  to  a  3a/3L*  "  «  type  of  stability  limit.  A  review  of  certain  dsta 
from  Ref.  38  seems  to  support  this  Interpretation  (Fig.  4.3).  In  this 
work,  tests  st  low  pressures  yielded  vslues  of  a  decreasing  to  zero  ss 
L*  Increased,  suggestive  of  curve  s  in  Fig.  4.1b.  As  pressure  wss 
raised,  it  became  increasingly  difficult  to  obtain  low  vslues  of  a, 
suggestive  of  s  trend  towards  s  stability  solution  like  curve  c  in 
Fig.  4.1b  (under  these  conditions,  when  testing  in  the  L*  rsnge  where 
a  -*■  0  on  the  lower  branch  of  curve  c,  oscillations  would  be  expected  to 
be  observed  lnstesd  st  the  higher  frequency  snd  a  corresponding  to  the 
upper  brsnch  of  the  a-L*  curve) .  The  behsvlor  at  higher  pressure  is 
thus  suggestive  of  stsblllty  limit  behavior  corresponding  to  3a/3L*  ■  «•, 
not  previously  considered  In  the  literature.  The  trend  from  sn  a  ■*  0 
to  a  3o/3L*  ■  «  type  of  stsblllty  limit  ss  pressure  Increases  esnnot  be 
encompassed  within  the  QSHOD  model-dependent  analyses  because  the  corre¬ 
sponding  form  of  the  A  versus  L*  curves  (Fig.  4.2)  are  relatively  insensi¬ 
tive  to  pressure  in  the  QSHOD  model. 
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FIG.  4.3.  Data  from  BMI  Teata  on  A-146  Propellant  (75X  Bimodal 
AP,  25 X  Polyurethane  Binder).  Trend  of  data  suggesta  transition 
from  a  «  0  type  atabllity  limit  for  low  frequency  (low  pressure , 
curve  A  (Fig.  4.1b))  teat* ,  to  3a/DL*  »  •  type  of  stability  limit 
in  high  frequency  (high  pressure t curve  B  (Fig.  4.1b))  tests.  Low 
values  of  a  could  not  be  obtained  at  high  frequency  or  pressure. 


FIG.  4.4.  Teat  Record  With  Transition  In  Oscilla¬ 
tion  Frequency  Suggestive  of  L*  Increasing  Through 
Crosa-over  Point  In  Stability  Curve  of  Type  C 
(Fig.  4.1b). 
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Figure  4.4  provides  an  interesting  (and  previously  enigmatic)  test 
of  the  perturbation  equations  (Eq  4.1  and  4.2).  .'.n  this  test,  as  the  L* 

increased  during  burning,  the  system  exhibited  a  spontaneous  abrupt  trans¬ 
ition  from  high  frequency  to  low  frequency  oscillations.  This  behavior 
is  easily  explained  in  terms  of  curve  b  of  Fig.  4.1b,  where  severe  high 
frequency  oscillations  are  to  be  expected  for  low  1*.  A  3a/3L*  •»  ■»  type 
"stability"  limit  is  reached,  beyond  which  oscillations  would  be  expected 
at  some  lower  frequency  corresponding  to  the  lower,  high  L*  section  of  the 
solution  curve.  From  an  examination  of  Fig.  4.2,  it  appears  that  this  be¬ 
havior  is  not  consistent  with  the  QSHOD  combustion  perturbation  models  be¬ 
cause  for  those  models  the  lower  branches  of  the  solution  curves  are  at 
negative  a.  It  is  the  "anomaly"  in  the  T  versus  u  curve  (in  Fig.  4.1a) 
that  causes  the  multiple  valued  functions  in  Fig.  4.1b.  In  view  of  the 
similar  qualitative  nature  of  the  a-L*  curves  in  Fig.  4.2,  and  in  view  of 
test  results  such  ss  illustrated  by  Fig.  4.4,  it  seems  entirely  plausible 
that  solutions  to  Eq.  4.1  and  4.2  with  real  response  functions  msy  have 
the  form  of  curve  b  in  Fig.  4.1b. 

From  the  foregoing  it  is  evident  that  there  is  much  to  be  learned 
by  considering  the  entire  solution  to  Eq.  4.1  and  4.2,  instead  of  just  the 
stability  limits.  Further,  it  seems  clear  that  the  QSHOD  models  for  per¬ 
turbation  combustion  are  qualitatively  incorrect.  Fecsuse  of  the  limita¬ 
tions  of  the  QSHOD  model,  attention  has  been  turned  to  a  more  detached 
study  of  the  structure  of  the  combustion  zone.  Motion  pictures  of  burning 
samples,  and  SEM  pictures  of  quenched  samples  have  revealed  in  consider¬ 
able  detail  the  nature  of  the  burning  surface  (see  earlier  section  of  this 
report).  It  consists  (at  low  pressures  typical  for  EMI)  of  an  array  of 
oxidizer  particles  protruding  from  the  pyrolyzing  binder.  It  must  be  pre¬ 
sumed  that  under  oscillatory  conditions,  these  oxidizer  particles  will  re¬ 
spond  in  a  different  way  than  will  the  binder,  with  resultant  oscillation 
in  mixture  ratio  of  gases  coming  from  the  combustion  zone.  Such  processes 
are  not  encompassed  in  the  QSHOD  perturbation  models  of  combustion.  Some 
success  (Ref.  38)  has  been  realized  in  an  alternate  scheme  for  correlating 
BM1  stability  limit  data,  a  scheme  in  which  the  natural  time  constant  for 
combustion  oscillations  is  assumed  to  be  the  time  to  repopulate  the  burning 
surface  with  oxidizer  particles  once  the  population  has  been  perturbed. 

This  time  constant  is  the  time  to  burn  the  distance  of  an  oxidizer  particle 
diameter.  For  some  propellants,  the  BM1  data  has  been  better  correlated  as 
a. function  of  pressure  and  oxidizer  particle  size  by  this  "layer"  concept 
ban  by  QSHOD  models. 

The  view  that  the  OSHOD  models  are  sometimes  even  qualitatively  in¬ 
accurate  is  supported  n  only  by  the  foregoing  examples,  but  also  by  sub¬ 
stantial  additional  experimental  evidence  (Ref.  11,  38,  46  and  47).  Con¬ 
tinued  reliance  upon  them  appears  to  be  blocking  full  and  correct  inter¬ 
pretations  of  BMI  experiments.  The  analytical  convenience  of  the  QSHOD 
model  is  unfortunately  also  an  intellectual  strait jacket.  It  may  be  that 
reproducibility  of  BMI  tests  is  too  poor  to  make  full  Interpretation  satis¬ 
fying,  or  the  applicability  of  Eq.  4.1  and  4.2  may  merit  some  scrutiny,  but 
the  most  promising  line  of  inquiry  seems  to  be  an  eacape  from  OSHOD  model- 
dependent  stability  limit  correlations  to  a  more  complete  and  less  inhibited 
interpretation  of  results. 
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5.  ACOUSTIC  INSTABILITY 


Studies  of  wave  mode  instability  have  continued  along  three  lines. 
The  dynamic  response  of  the  combustion  zone  to  peipendicularly  incident 
pressure  waves  has  been  studied  by  both  analytical  modeling  and  measure¬ 
ments  of  detailed  combustion  zone  behavior.  The  net  effect  of  this 
press ure-cou,  led  combustion  response  in  amplifying  pressure  waves  was 
studied  experimentally  using  the  T-burner,  with  several  propellent  formu¬ 
lations.  Finally,  a  continuing  effort  wss  devoted  to  testing  the 
pressure-coupled  response  of  propellants  used  in  development  programs, 
with  the  purposes  of  providing  advance  warning  of  risk,  and  evaluation 
of  relevance  of  the  laboratory  tests  to  observed  motor  behavior. 


5.1.  ANALYTICAL  MODELING  OF  PRESSURE  COUPLING 

A  survey  of  one-dimensional  models  of  pressure-coupling  was  pre¬ 
sented  in  an  earlier  report.  This  work  was  rewritten,  updated,1  and  to 
some  extent  simplified  during  the  present  reporting  period,  snd  was 
published  as  a  survey  paper  in  a  technicsl  journal  (Ref.  39).  This  work 
was  more  successful  than  originally  anticipated,  .n  thst  it  provided  the 
principal  -esults  of  all  previous  models  in  terms  of  a  single  functional 
form  for  the  burning  rate  response  function. 


m 1  /m 
p’/p 


n  AB  +  n  (X-l) 

_  S 

X  +  (A/X)  -  (1  +  A)  +  AB 


(5.1) 


In  this  expression  the  left-hand  side  is  essentislly  the  amplification 
factor  for  incident  pressure  waves.  The  values  of  the  quantities  A,  B 
end  ns  depend  upon  fundamental  aspects  of  the  combustion  such  as  acti¬ 
vation  energies,  and  X  is  s  function  of  the  normalized  frequency  of  os¬ 
cillation.  From  this  relation  it  is  easy  to  determine  the  relative  con¬ 
tribution  of  the  various  reaction  steps  to  instability  in  terms  of  the 

psrameters  A,  B,  and  ng,  assuming  those  parameters  are  known  More 

recent  work,  described  in  this  report,  has  been  focused  on  elucidation 
of  the  processes  determining  the  parameters  A,  B,  and  ns  and  on  develop¬ 
ment  of  more  rigorous  modeling  of  the  combustion  response. 


1  Work  funded  psrtly  from  Naval  Ordnance  Systems  Command  Task  No. 
UF  19.332.330. 
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One  of  the  more  conspicuous  aspects  of  the  models  represented  by 
Eq.  5.1  ia  a  "persistent"  prediction  that  the  maximum  value  of  the  re¬ 
sponse  function  (m'/m)/(p'/p)  la  relatively  Insensitive  to  pressure, 
and  that  the  maximum  occurs  at  a  frequency  of  preseure  oscillation  that 
is  a  function  of  burning  rate, 


This  point  is  mentioned  here  because  the  experiaentil  results  reported 
below  are  in  sharp  contrast  to  this  prediction. 

While  it  wilj  no  doubt  be  necessary  to  go  to  more  complicated 
modeling  of  the  combustion  rone  if  accurate  predictive  capability  is  to 
be  achieved,  there  are  some  rather  elementary  extensions  sf  the  present 
models  that  would  enhance  their  relevance.  One  of  tr.a  principal  weak¬ 
nesses  of  previous  models  has  been  in  representation  of  the  beet  trans¬ 
fer  to  the  burning  surface  from  the  gas  phase.  This  has  in  most  cases 
been  contained  in  the  undetermined  parameter  "B"  in  Eq.  5.1,  and  in  all 
cases  has  been  based  on  the  assumption  that  the  response  of  the  gas 
phase  during  oscillations  was  quasi-static.  In  one  case  where  B  was 
not  treated  in  a  completely  parametric  way,  it  was  assumed  (Ref.  11) 
that  the  gas  phase  reaction  started  at  the  surface,  and  remained  at 
constant  rate  until  the  reactants  were  consumed.  The  analytical  model 
was  extended  during  this  reporting  period  to  encompass  s  more  realistic 
representation  of  the  distribution  of  the  gas  phase  reactions,  and  a 
calculation  was  also  made  of  the  limitations  of  the  quasi-static  approx¬ 
imation  for  the  gas  phsse  combustion.  These  developments  are  reported 
fully  in  Ref.  18.  The  results  indicate  that  the  surface-attached  gas 
phase  reaction  zone  model  used  by  Krier,  et  al,  (Ref.  52)  gives  signifi¬ 
cantly  different  results  than  a  model  in  which  some  stand-off  of  the 
flame  is  allowed.  Examinstion  of  the  effect  of  the  quasi-static  assump¬ 
tion  shows  that  it  does  not  materially  affect  the  prediction  of  the  in- 
phase  component  of  combustion  response  to  pressure  in  the  frequency 
range  up  to  10,000  cps,  but  does  cause  error  in  calculation  of  tbe  magni¬ 
tude  and  phase  of  the  overall  combustion  response  at  high  frequency. 
These  results  presumably  improve  the  prospects  of  correlation  of  experi¬ 
ments  with  theory,  but  the  experimental  work  on  combustion  zone  struc¬ 
ture  indicates  that  still  more  complicated  models  will  be  required  be¬ 
fore  fully  satisfactory  agreement  will  be  obtained  between  experiment 
and  theory.  Existing  models  may  have  some  predictive  capability  for 
trends.  One  remaining  extension  of  the  one-dimensions 1  theory  is  the 
nonlinear  analysis,  which  is  currently  in  progress. 

5.2.  MAPPING  RESPONSE  FUNCTION  OF  PROPELLANTS 

One  of  the  long-time  gosls  of  experimental  research  on  combustion 
instability  has  been  to  characterize  the  perturbation  behavior  of  pro¬ 
pellants  over  the  entire  domain  of  pressure  and  frequency  for  which 


ic  'ability  1$  possible.  In  the  preseat  program,  extensive  effort  has 
been  devoted  to  development  of  experimental  techniques  for  this  purpose, 
but  eacb  technique  has  been  amenable  to  only  a  restricted  range  of  con¬ 
ditions  of  interest.  During  the  past  year,  tests  have  been  completed 
on  certain  propellants  over  the  range  of  conditions  accessible  to  three 
experimental  systems.  The  response  function  was  measured  in  the  bigh- 
frequency  domain  using  the  1.5-inch  diameter  T-burner  (Ref.  47  and  49), 
nnd  in  tne  low  frequency  domain  using  the  5.5-lnch  diameter  T-bumer 
(Ref.'  1,  2,  4,  5,  50  and  51),  and  the  l.^-inch  diameter  L*  burner  (Ref. 
1-5,  17,  38  and  50).  The  results  of  these  studies  were  reported  at  the 
Twelfth  Combustion  Symposium  (Ref.  46),  and  will  appesr  in  the  published 
volume  of  the  meeting.  The  results  are  summarized  briefly  here. 

Studies  were  made  of  four  propellants,  designated  A-13,  A-35,  A-91 
awl  A-148.  A-13  consists  of  a  76-242  mixture  of  80  v  AP  and  PBAN  binder, 

A-35  consists  of  a  75-252  mixture  of  80  y  AP  and  polyurethane  binder, 
and  A-91  consists  of  a  67-25-82  mixture  of  80  w  AP,  polyurethane  binder 
and  80  p  aluminum.  A-148  is  a  75-252  mixture  of  bimodal  AP  (50/50  mix¬ 
ture  of  15  p  end  200  p  particle  sizes)  and  polyurethane  binder.  The 
test  results  were  summarized  in  the  form  of  contour  maps  of, the  in-phase 
component  of  the  response  function  (p/e),  plotted  in  the ■  pressure- 
frequency  plane. 

Recalling  the  discussion  in  the  last  section,  the  one-dimensional 
perturbation  theories  predict  that  the  response  function  ''surface" 
above  the  p-f  plane  should  be  a  ridge  of  roughly  constant  height 
(Fig.  5.1a),  with  the  corresponding  constant-y/e  lines  in  the  p-f  plane 
being  open  curves  following  the  trend  of  the  ridge  of  the  response  func¬ 
tion  surface.  The  experimental  results  did  not  conform  to  this  trend 
(Fig.  5.2),  indicating  the  need  for  more  realistic  modeling  of  the  com¬ 
bustion  process.  All  four  propellants  showed  behavior  suggestive  of  a 
"peaked"  response  function  surface  (Fig.  5.1b)  instead  of  a  ridge,  with 
the  location  and  height  of  the  peak  depending  on  the  propellant  com¬ 
position.  A-13  propellant  (PBAN  binder)  showed  a  broad  peak  centered 
around  300  psia  and  300  cps  (Fig.  5.2a).  A-35  propellant  (polyurethane 

binder)  showed  a  peak  centered  around  70  psia  and  60  cp3  (only  roughly 
determined),  with  magnitude  comparable  to  that  of  Fig.  5.2b.  The 
"bimodal"  formulation  A-148  gave  similar  trends  (Fig.  5.2c)  with  no 
singular  behavior  currently  attributable  to  the  bimodal  oxidizer  parti¬ 
cle  size  distribution.  The  response  function  contours  are  open  in 
Fig.  5.2c  because  the  unstable  region  extended  into  the  region  of 
atmospheric  pressure  where  testing  was  not  feasible.  The  polyurethane 
formulation  with  82  aluminum  (A-91)  showed  a  very  high,  localized  peak 
in  the  region  of  30  psia  and  15  ops  (Fig.  5. 2d) ,  a  result  that  con¬ 
firmed  previous  qualitative  interpretations  to  the  effect  that  the  metal 
combustion  can  produce  a  "preferred  frequency"  type  of  instability 
(Ref.  1-4  and  51)  that  can  be  quite  severe  under  rather  specialized  con¬ 
ditions.  The  results  also  indicate  that  instability  problems  tend  to 
be  more  serio  at  low  frequencies  and  low  pressures,  although  results 
with  other  pi  .  Hants  indicate  that  the  examples  shown  here  "peak  out" 
at  lower  pres  re  and  frequency  than  the  average. 


FIG.  5.1.  (a)  Response  Function  Surface  Above  the  Pressure- 
Frequency  Plane  as  Predicted  by  One-dimensional  Perturbation 
Theories,  (b)  Idealized  Representation  of  the  Response  Func 
tion  Surface  Abo"e  the  Pressure-Frequency  Plane  as  Determined 
by  Experimental  Measurement. 
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5.3.  OBSERVATIONS  OF  COMtJSTION  ZONE  DYNAMICS 


Direct  observation  of  combustion  zone  dynamics  is  seldom  possible, 
yet  knowledge  of  combustion  zone  dynamics  is  vital  to  an  understanding 
of  phenomena  involved  in  acoustic  instability  and  to  the  construction 
of  correct  mathematical  models.  Most  of  our  knowledge  of  combustion 
zone  dynamics  is  obtained  through  indirect  evidence.  The  NVC  5.5-inch 
diameter  DECEV  (T-burner)  ia  being  used  as  a  tool  to  delve  into  certain 
aspects  of  combustion  zone  dynamics.  Although  events  occurring  in  the 
combustion  zone  are  Inferred  from  observations  of  the  combustion  products 
as  they  flow  away  from  the  propellant  surface,  several  sdvances  have  been 
made  in  understanding  these  complicated  phenomena.  The  5.5-inch  DECEV 
haa  been  Instrumented  with  a  variety  of  devices  in  addition  to  the  usual 
pressure  transducers,  which  Include  multiple  radiation  transducers  to 
detect  metal  combustion  clouds  (Ref.  4),  a  slit  window  and  camera  ror 
space-time  recording  of  gaa  motion  (Ref.  5) ,  and  high  speed  gas  tempera¬ 
ture  measurements  (Ref.  5). 

Earlier  difficulties  were  encountered  in  interpreting  streak  pho¬ 
tographs  of  the  gas  motion  in  the  burner  (Ref.  5).  These  were  largely 
surmounted  during  the  present  reporting  period  by  moving  the  propellant 
to  a  position  initially  at  one  end  of  the  slit  window,  by  correcting  the 
image  exposure  on  the  film,  and  by  completing  the  development  of  a  mathe¬ 
matical  model  which  describes  gss  flow  in  the  burner  (Ref.  6).  The 
mathematical  model  was  programmed  for  digital  computer  and  solutions  com¬ 
pared  with  streak  photographs  of  gas  motion  in  the  burner.  The  compari¬ 
son  shows  that  the  model  provides  a  reasonably  accurate  description 
(see  Fig.  4.1  in  Ref.  6). 

Early  streak  data  indicated  that  temperature  variations  in  the  flow 
arose  from  an  apparent  flame  temperature  variation  as  well  as  the  pre¬ 
dicted  entropy  variations  (Ref.  5).  Measurements  of  temperature  oscil¬ 
lations  near  the  propellant  surface  gave  results  that  in  some  cases  were 
much  larger  than  could  be  accounted  for  by  isentroplc  compression  of  the 
gas.  Although  the  grester  values  might  be  explained  by  a  dependence  of 
the  flame  temperature  on  pressure,  measurements  of  the  temperature  under 
nonoscillating  conditions  anH  calculation  of  the  dependence  of  temper¬ 
ature  from  equilibrium  data  tend  to  discount  this  explanation.  Another 
explanation  was  sought  along  the  line  of  changing  gas  stoichiometry 
caused  by  differential  changes  in  burning  rate  of  the  propellant  con¬ 
stituents  with  pressure.  A  series  of  calculations  was  made  on  a  digital 
computer  for  a  hypothetical  propellant  with  several  fuel  to  oxidizer 
ratios  and  chamber  pressures.  The  calculated  flame  temperatures  are 
plotted  in  Fig.  5.3.  Several  significant  conclusions  can  be  drawn  from 
the  data:  (1)  The  temperature  is  most  sensitive  to  the  fuel-oxidizer 
ratio  in  the  region  representing  propellants  being  investigated  (75* 
oxidizer),  (2)  the  temperature  is  least  sensitive  to  changes  in  mean 
pressure  in  this  same  region,  (3)  the  sensitivity  of  temperature  to 
changes  in  fuel-to-cxidlzer  ratio  in  the  70*  to  80*  oxidizer  region  is 
sufficient  (over  100°X  for  a  one  percent  oxidizer  change)  to  explain 
the  excessive  temperature  excursions  noted  earlier. 
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FIG,  5.3.  Dependence  of  Flame  Temper¬ 
ature  on  Percent  Oxidizer  and  on  Pres- 
sure  for  NH^CIO^-Estane  Propellants. 


Numerous  problems  remain  unsolved  in  this  subject  area.  Some  of 
the  points  mentioned  below  are  Interrelated  but  a  partial  list  Includes: 

1.  Attempt  characterization  of  the  combustion  zone,  particularly 
its  extent  into  the  flow  field,  which  may  be  considerable  when  dealing 
with  metallized  propellants. 

2.  Measure  temperatures  at  known  distances  from  the  propellant  sur¬ 
face  to:  (a)  verify  that  the  theoretical  model  of  gas  flow  in  the  burner 
provides  an  adequate  description  of  events;  (b)  use  data,  when  point  (a) 
is  confirmed,  to  extrapolate  back  to  the  combustion  zone  to  calculate 
flame  behavior. 

3.  Gather  evidence  to  support  or  deny  "mixture  ratio"  hypothesis. 

4.  Apply  standard  and  time-resolved  spectroscopy  to  study  optical 
emissions  from  combustion  gases  to:  (a)  identify  general  nature  of 
emission,  (b)  study  the  nature  of  reactions  in  the  comoustion  zone,  (c) 
determine  the  spatial  extent  of  combustion  zone,  (d)  map  distribution  of 
chemical  species  in  space  and  time  and  attempt  to  infer  combustion  zone 
dynamic  behavior. 
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5.  Extend  atudiea  to  higher  preasurea  and  frequencies. 

6.  Develop  more  adequate  means  for  study  of  metal  combust!  in  and 
the  role  of  metal  In  combuation  oscillatlona. 


5.4.  TESTING  OF  PROPELLANTS  FROM  DEVELOPMENT  PROGRAMS 

Four  propellants  from  the  Jet  Propulsion  Laboratory  were  evaluated 
for  combustion  characteristics.  These  prnpellants  contained  a  "sature- 
tbane"  binder  and  were  designed  to  withstand  the  high  temperature  con¬ 
ditions  necessary  to  biologically  sterilize  the  apparatus  containing  the 
pnpellant .  They  were  evaluated  in  the  photographic  strand  bomb  and  in 
the  NWC  1.5-inch  diameter  T-bumer.  Two  of  the  compositions  were  non- 
meta.'lized,  one  contained  1/2Z  powdered  aluminum,  and  the  fourth  had  2Z 
alumj  num. 

Color  photography  of  the  combustion  zone  showed  somewhat  mor  molten 
binder  than  with  NWC  polyurethane  compositions.  T-burner  testing  on  the 
propellants  showed  moderate  instability  for  the  two  unmetallized  com¬ 
positions  in  the  mid-range  of  pressure  and  frequency  (see  Fig.  5.4). 
Acoustic  response  function  values  were  moderately  low  in  the  range  teated 
and  acoustic  instability  appears  unlikely  in  the  motor  design  currently 
being  considered  by  this  propellant. 

A  few  preliminary  tests  in  the  L*-burner  resulted  in  oscillatory 
behavior  at  extremely  low  L*  and  pressure.  Nonacoustic  instability 
would  seem  unlikely  in  the  L*-pressure  range  of  the  current  motor  design. 

Two  compositions  were  received  from  Thiokol  Chemical  Co.  and  were 
tested  in  the  NWC  1.5-inch  T-burner.  Both  showed  moderate  instability. 
One  of  these  was  a  heavily  metallized  composition  with  a  high-energy 
oxidizer.  In  the  four  cases  where  this  propellant  was  unstable,  the 
response  function  was  of  the  order  of  unity  or  less  in  three  tests  and 
was  less  than  two  on  the  other  test.  These  are  only  moderate  values  of 
the  response  function,  in  spite  of  the  fact  that  the  growth  and  decay 
rates  of  the  oscillations  were  of  somewhat  greater  magnitude  than  usu¬ 
ally  encountered.  A  relatively  high  flame  temperature  and  burning  rate 
appear  to  be  the  principal  factors  causing  the  instability  of  this  pro¬ 
pellant,  which  is  capable  of  producing  high  growth  a's  even  in  the 
presence  of  the  high  attenuation  (high  decay  a)  resulting  from  the 
aluminum  oxide  smoke.  The  likelihood  of  acoustic  instability  with  this 
propellant  in  the  proposed  application  was  judged  to  be  moderately  low 
for  axial  modes  and  somewhat  higher  for  trans/erse  modes.  As  in  all 
cases,  precise  prediction  of  motor  stability  from  T-burner  results  is 
not  possible  because  the  gain-loss  characteristics  of  the  motor  are  ex¬ 
tremely  difficult  to  determine  and  the  influence  of  gas  perturbations 
parallel  tn  the  burning  surface  is  not  tested  in  the  T-burner. 
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FIG.  5.4-  Acoustic  Instability  Behavior  of  a  JPL  "Saturethane" 
Propellant.  Open  symbols  represent  stable  tests,  solid  symbols 
denote  unstable  behavior.  Numbers  indicate  values  of  the  in- 
phase  component  of  the  response  function.  The  approximate  re¬ 
gion  of  T-burner  lnatabillty  with  this  propellant  is  Indicated 
by  the  shaded  area. 


As  a  result  of  instability  troubles  in  a  motor  under  development, 
six  candidate  propellant  compositions  were  sent  by  Aerojet-General  Corp., 
which  were  tested  in  the  5.5-inch  diameter  T-burner  at  a  nominal  fre¬ 
quency  of  380  cps.  Only  three  tests  of  each  propellant  were  made  due 
to  shortage  of  material.  Two  tests  series  were  at  ambient  temperature 
and  the  final  one  was  with  propellant  at  approximately  150°F.  No  insta¬ 
bility  was  noted  for  any  of  the  compositions  at  any  condition.  Pho¬ 
tographic  strand  bomb  resulta,  however,  showed  a  wide  variety  of  alu¬ 
minum  behavior.  It  seems  likely  that  the  instability  encountered  in 
the  motor  was  due  primarily  to  velocity-coupled  mechanisms  rather  than 
to  pressure-coupled  phenomene. 
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6.  ACOUSTIC  ATTENUATION 


Stability  of  solid  propellant  combustors  is  primarily  s  consequence 
of  attenuation  of  flow  disturbances  by  viscous  processes  snd  by  convec¬ 
tion  through  the  nozzle.  Because  of  the  difficulty  of  calculating  these 
processes  in  combustors  of  complicated  geomet.y,  a  cold-flow  experimen¬ 
tal  technique  was  developed  which  would  permit  determination  of  the 
effect  of  geometrical  variables  on  losses  in  axial  modes  of  combustor 
models.  This  work  was  reported  previously  (Ref.  8,  9,  and  53).  Experi¬ 
mental  work  during  the  present  reporting  period  has  been  limited  to 
attenustion  measurements  incidental  to  T-burner  testing,  the  cold-flow 
work  having  been  discontinued  due  to  reduced  funding. 

An  analytical  study  was  completed  concerning  the  interpretation  of 
cold-flow  experiments.  This  study,  entitled  "An  Analysis  of  Axial 
Acoustic  Waves  in  s  Cold-Flow  Rocket",  hss  been  reported  fully  in  Ref. 

10  and  will  be  discussed  only  briefly  here.  In  the  study,  equations 
were  derived  for  steady  state  or  driven  oscillations  and  for  freely  de¬ 
caying  waves.  The  analysis  was  motivated  by  an  uncertainty  regarding 
the  quantitative  validity  of  the  "half  power  bandwidth"  method  used  in 
determining  the  attenuation  alpha  from  resonance  curves  obtained  from 
the  cold-flow  tests.  The  analysis  showed  that  the  method  was  applicable 
to  cases  where  the  effective  admittance  of  the  nozzle  end  of  the  motor 
was  small  snd  the  motor  geometry  relatively  simple.  For  more  "lossy" 
motors  an  alternative  method  was  provided  in  which  the  acoustic  attenu¬ 
ation  may  be  obtained  graphically  from  plots  of  the  effective  nozzle 
admittance  dependence  on  resonant  frequency  and  half  power  bandwidth. 
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7.  CONSULTING  AND  REPORTING 


An  appreciable  amount  of  time  was  apent  during  this  reporting 
period  on  various  consulting  and  reporting  activities  that  would  be 
regarded  as  "bridging  the  gap"  between  research  and  application.  These 
a- livltles  are  summarized  briefly  below: 

7.1.  CONSULTING  ON  DEVELOPMENT  PROGRAMS 

A  review  of  development  and  proposed  application  of  aaturethane 
propellants  was  made  for  the  Jet  Propulsion  Laboratory/California  Insti¬ 
tute  of  Technology  (JPL/CIT),  relative  to  the  llkllhood  of  problems  with 
combuation  instability.  Tests  were  made  of  samples  of  propellant 
supplied  by  JPL/CIT,  tests  including  window  bomb  photography,  T-burner 
firings  and  nonacoustic  burner  tests.  From  the  results  of  these  studies 
it  was  concluded  that  the  susceptibility  of  the  saturethane  propellants 
tested  to  combustion  is  relatively  low;  nonacoustic  instability  Is  a 
possibility,  but  only  at  lower  pressure  and  L*  than  called  for  In  preaent 
designs . 

Oscillatory  behavior  (axial  mode)  In  developmental  static  firings 
of  the  motor  for  the  Navy  Standard  Missile  was  reviewed  and  recommenda¬ 
tions  made  for  modification  of  oxidizer  and  aluminum  particle  size  In 
the  propellant  and  a  minor  change  In  charge  design.  The  changes  led  to 
stable  behavior. 

Axial  mode  oscillatory  behavior  in  a  Rockrtdyne  advanced  motor  for 
Redeye  was  studied  and  burner  tests  were  conducted  on  a  family  of  pro¬ 
pellants  with  variations  In  oxidizer  and  aluminum  particle  alze  dis¬ 
tribution.  The  objective  of  these  studies  (Insofar  as  the  NASA  p-'-gram 
is  concerned)  was  to  determine  the  relevance  of  laboratory  tests  to 
actual  motor  behavior.  This  particular  motor  represents  n  motor  test 
situation  (aluminized  propellant,  axial  mode  oscillation,  frequency 
600  cps)  In  which  only  limited  success  has  been  echleveJ  in  laboratory 
testing  methods.  In  this  case,  the  window  bomb  photography  showed  con¬ 
siderable  difference  in  steady-state  combustlcn  behavior  among  the  pro¬ 
pellants  tested,  but  motor  testing  was  confined  (for  separate  reasons) 
to  a  subaet  of  propellant  formulations  whose  combustion  in  the  window 
bomb  appeared  similar  Al’  were  unstable  In  rocket  motor  firings.  Ex¬ 
tensive  testing  Is  being  done  In  a  5.5-lnch  diameter  T-burner.  Oscil¬ 
lations  do  not  occur  spontaneously,  presumably  because  cf  high  attenu¬ 
ation  by  AI2O3  smoke.  A  pulse-decay  technique  used  previously  for  low 
frequency  Instability  testing  Is  being  adapted  to  this  600  cps  range. 
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In  cooperation  with  Thiokol  Chemical  Corporation  and  the  Air  Porce 
Rocket  Propulsion  Laboratory,  tests  (T-turner)  were  run  on  several  pro¬ 
pellant  formulations  using  hydrazine  dlperchlorste  as  the  oxidizer.  The 
testing  work  was  done  on  s  Navy  program,  but  Is  mentioned  here  because 
of  its  fundamental  Importance.  These  oropellsnts  are  unusually  unstable 
at  the  high  frequencies  tested,  and  formulations  with  up  to  4Z  aluminum 
were  still  very  unstable.  Although  this  and  other  test  work  hss  shown 
a  serious  problem  with  Instability,  it  should  be  noted  that  no  informa¬ 
tion  has  yet  been  obtained  at  the  frequencies  characteristic  of  motors 
of  operational  size.  Thus  the  uufavorable  results  may  lead  to  unjusti¬ 
fied  concern  with  this  problem.  In  another  respect  the  combustion  insta¬ 
bility  problem  with  HP2  is  of  continuing  interest:  why  are  HP2  propel¬ 
lants  so  unstable?  The  answer  to  that  question  may  reveal  much  sbout 
the  detailed  mechanisms  of  Instability. 

Another  development  program  that  has  encountered  s  significant  de¬ 
gree  of  oscillatory  combustion  Is  the  Mlnuteman  3rd  Stage  Wing  XI,  This 
experience  was  reviewed  with  personnel  of  Hercules  Powder  Company,  and 
discussions  are  continuing.  The  behavior  Involves  axlsl  mode  oscilla¬ 
tions  at  about  500  cps.  The  principal  interest  to  NASA  is  the  fact  that 
this  motor  Is  similar  to  high  performance  "space"  motors  NASA  might  use, 
with  s  heavily  aluminized  propellant,  lov  operating  pressure,  and  high 
volumetric  loading. 

Axial  mode  oscillatory  behavior  was  encountered  In  static  develop¬ 
ment  tests  on  an  advanced  motor  for  Sparrow,  st  frequencies  in  t'-.e  400- 
1200  cps  range.  In  consultations  on  this  program  it  was  recommended 
that  s  Helmholtz  resonator  attenuator  be  used  in  the  motor.  Tests  to 
date  have  shown  the  attenuator  to  be  effective. 

A  general  comment  on  the  above  and  OLher  recent  experience  with 
combustion  Instability  Is  needed.  Mopt  recent  problems  have  been  In  the 
100-1200  cps  range.  In  axlsl  modes,  with  heavily  aluminized  propellants. 
There  Is  a  vide-privalling  misconception  that  large  motors  and  motors 
with  heavily  aluminized  propellants  are  Inherently  stable.  While  most 
Information  supports  these  views  relative  to  the  big  booster  motors.  It 
Is  absolutely  clear  that  It  Is  not  true  for  intermediate  size  boosters 
and  upper  stages.  There  are  twe  programs  in  which  oscillatory  behavior 
at  60  cps  hss  occurred.  In  other  development  programs  It  has  been  re¬ 
peatedly  observed  that,  this  sxlal  mode  Instability  can  occur  consistently 
in  rather  Innocuous  form  and  then  undergo  transition  to  violently  un¬ 
stable  behavior  because  of  small  propellant  or  design  changes,  or  flow 
disturbances  due  to  discharge  of  debris  through  the  nozzle.  Thus  in¬ 
creasing  attention  to  this  class  of  behavior  is  needed. 

7.2.  SURVEY  AND  EDUCATIONAL  ACTIVITIES 

Efforts  have  continued  to  organize  knowledge  in  the  field  of  com¬ 
bustion  instability  and  present  it  in  s  form  conducive  to  more  widespread 
application.  A  survey  paper  was  prepared  and  presented  at  the  Twelfth 
Combustion  Symposium  (Ref.  54),  summarizing  recent  developments  in 
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research.  A  tutorial  paper  was  preiented  at  the  Third  ICRPC/AIAA  Solid 
Propulsion  Meeting,  8 unsnarl ing  the  ncture,  prevalence  and  scientific 
basis  of  combustion  instability,  and  this  paper  le  being  expanded  for 
publication  through  CFIA.  A  two-d*1*  seminar  was  presented  on  the  topic 
of  solid  propellant  combustion,  under  the  ausplcee  of  the  American  Instl- 
tutute  of  Aerouautlca  snd  Astronautics.  Thie  couree  drew  heavilv  on 
results  from  the  preeent  program,  including  e  tutorial  preeentetlon  of 
ingredient  decompoeltion,  et early  etate  combustion,  instability  end  igni¬ 
tion.  Coplea  of  the  courae  notea  are  available  from  AIAA  Headquarters 
(Ref.  28).  The  material  is  currently  being  rewritten  snd  expended  into 
a  book. 

7.3.  SUPPORT  OF  THE  INTERAGENCY  CHEMICAL 
ROCKET  PROPULSION  GROUP 

Activities  in  this  erea  have  Included  chairmanship  of  the  Solid 
Propellant  Combustion  Working  Group,  organlzetion  of  a  tutorial  session 
at  the  Third  ICRPG/AIAA  Solid  Propulsion  Meeting,  membership  in  snd 
eupport  of  the  "T-Buroer  Standardization"  Committee  snd  the  "Arc-Image 
Ignition  Standardization"  Committee,  end  Initial  planning  of  a  working 
meeting  on  ammonium  perchlorate  decomposition  (for  April  30,  1969). 

7.4.  SUMMARY  OP  REPORTS 

A  list  of  publications  pertaining  to  work  on  tbis  program  is 
enclosed  as  Appendix  A. 
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